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Abstract 

The entanglement entropy is a fundamental quantity which characterizes the correlations 
between sub-systems in a larger quantum-mechanical system. For two sub-systems separated 
by a surface the entanglement entropy is proportional to the area of the surface and depends on 
the UV cutoff which regulates the short-distance correlations. The geometrical nature of the 
entanglement entropy calculation is particularly intriguing when applied to black holes when 
the entangling surface is the black hole horizon. 1 review a variety of aspects of this calculation: 
the useful mathematical tools such as the geometry of spaces with conical singularities and 
the heat kernel method, the UV divergences in the entropy and their renormalization, the 
logarithmic terms in the entanglement entropy in 4 and 6 dimensions and their relation to the 
conformal anomalies. The focus in the review is on the systematic use of the conical singularity 
method. The relations to other known approaches such as 't Hooft's brick wall model and the 
Euclidean path integral in the optical metric are discussed in detail. The puzzling behavior of 
the entanglement entropy due to fields which non-minimally couple to gravity is emphasized. 
The holographic description of the entanglement entropy of the black hole horizon is illustrated 
on the two- and four-dimensional examples. Finally, I examine the possibility to interpret the 
Bekenstein-Hawking entropy entirely as the entanglement entropy. 
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1 Introduction 



One of the mysteries in modern physics is why black holes have an entropy. This entropy, known 
as the Bekenstein-Hawking entropy, was first introduced by Bekenstein ,3T], [T5], [TH] rather as a 
useful analogy. Soon after that, this idea was put on a firm ground by Hawking |129| who showed 
that black holes thermally radiate and calculated the black hole temperature. The main feature of 
the Bekenstein-Hawking entropy is its proportionality to the area of the black hole horizon. This 
property makes it rather different from the usual entropy, for example the entropy of a thermal 
gas in a box, which is proportional to the volume. 

In 1986 Bombelli, Koul, Lee and Sorkin [53] published a paper where they considered the so- 
called reduced density matrix, obtained by tracing over the degrees of freedom of a quantum field 
that are inside the horizon. This procedure appears to be very natural for black holes, since the 
black hole horizon plays the role of a causal boundary, which does not allow anyone outside the 
black hole to have access to the events which take place inside the horizon. Another attempt to 
understand the entropy of black holes was made by 't Hooft in 1985 |213j . His idea was to calculate 
the entropy of the thermal gas of Hawking particles which propagate just outside the horizon. This 
calculation has uncovered two remarkable features: the entropy does turn out to be proportional 
to the horizon area, however, in order to regularize the density of states very close to the horizon, 
it was necessary to introduce the so-called "brick wall", a boundary which is placed at a small 
distance from the actual horizon. This small distance plays the role of a regulator in the 't Hooft's 
calculation. Thus, the first indications that entropy may grow as area were found. 

An important step in the development of these ideas was made in 1993 when a paper of 
Srednicki |207] appeared. In this very inspiring paper Srednicki calculated the reduced density 
and the corresponding entropy directly in flat spacetime by tracing over the degrees of freedom 
residing inside an imaginary surface. The entropy defined in this calculation has became known as 
the entanglement entropy. Sometimes the term geometric entropy is used as well. The entanglement 
entropy, as was shown by Srednicki, is proportional to the area of the entangling surface. This 
fact is naturally explained by observing that the entanglement entropy is non- vanishing due to the 
short-distance correlations present in the system. Thus only modes which are located in a small 
region close to the surface contribute to the entropy. By virtue of this fact one finds that the size 
of this region plays the role of the UV regulator so that the entanglement entropy is a UV sensitive 
quantity. A surprising feature of Srednicki's calculation is that no black hole is actually needed: 
the entanglement entropy of a quantum field in fiat spacetime already establishes the area law. In 
an independent paper Frolov and Novikov |99j applied a similar approach directly to a black hole. 
These results have sparked the interest in the entanglement entropy. In particular, it was realized 
that the " brick wall" model of t'Hooft studies a similar entropy and that the two entropies are 
in fact related. On the technical side of the problem, a very efficient method was developed to 
calculate the entanglement entropy. This method, first considered by Susskind [210] . is based on 
a simple replica trick in which one first introduces a small conical singularity at the entangling 
surface, evaluates the effective action of a quantum field on the background of metric with a conical 
singularity and then differentiates the action with respect to the deficit angle. By means of this 
method one has developed a systematic calculation of the UV divergent terms in the geometric 
entropy of black holes, revealing the covariant structure of the divergences [33], |195j . [112) . In 
particular, the logarithmic UV divergent terms in the entropy were found [197] . The other aspect, 
which was widely discussed in the literature, is whether the UV divergence in the entanglement 
entropy could be properly renormalized. It was suggested by Susskind and Uglum [212] that the 
standard renormalization of Newton's constant makes the entropy finite provided one considers 
the entanglement entropy as a quantum contribution to the Bekenstein-Hawking entropy. This 
proposal however did not answer the question whether the Bekenstein-Hawking entropy itself can 
be considered as an entropy of entanglement. It was proposed by Jacobson [142] that, in models 
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in which Newton's constant is induced in the spirit of Sakharov's ideas, the Bckenstein-Hawking 
entropy would also be properly induced. A concrete model to test this idea was considered in 

Unfortunately, in the 90-s the study of the entanglement entropy could not compete with the 
booming success of the string theory (based on D-branes) calculations of the black hole entropy 
[208]. The second wave of interest in the entanglement entropy has started in 2003 with works 
studying the entropy in condensed matter systems and in lattice models. These studies revealed 
the universality of the approach based on the replica trick and the efhciency of the conformal 
symmetry to compute the entropy in two dimensions. The black holes again became in the focus 
of study in 2006 after work of Ryu and Takayanagi [188] where a holographic interpretation of the 
entanglement entropy was proposed. In this proposal, in the frame of the AdS/CFT correspon- 
dence, the entanglement entropy, defined on a boundary of anti-de Sitter, is related to the area 
of a certain minimal surface in the bulk of the anti-de Sitter spacetime. This proposal opened 
interesting possibilities for computing, in a purely geometrical way, the entropy and for addressing 
in a new setting the question of the statistical interpretation of the Bekenstein- Hawking entropy. 

The progress made in the recent years and the intensity of the on-going research indicate that 
the entanglement entropy is a very promising direction which in the coming years may lead to a 
breakthrough in our understanding of the black holes and Quantum Gravity. A number of very 
nice reviews appeared in the recent years that address the role of the entanglement entropy for 
black holes ^U\, ^U\, |146J . [53] : review the calculation of the entanglement entropy in quantum 
field theory in flat spacetime [81], [37] and the role of the conformal symmetry [29]; focus on the 
holographic aspects of the entanglement entropy jl84J . [12j . In the present review I build on these 
works and focus on the study of the entanglement entropy as applied to black holes. The goal 
of this review is to collect a complete variety of results and present them in a systematic and 
self-consistent way without neglecting neither technical nor principal aspects of the problem. 



2 Entanglement entropy in Minkowski space-time 
2.1 Definition 

Consider a pure vacuum state > of a quantum system defined inside a space-like region O 
and suppose that the degrees of freedom in the system can be considered as located inside certain 
sub- regions of C A simple example of this sort is a system of coupled oscillators placed in the sites 
of a space-like lattice. Then for an arbitrary imaginary surface E which separates the region O in 
two complementary sub-regions A and -B, the system in question can be represented as a union of 
two sub-systems. The wave function of the global system is given by a linear combination of the 
product of quantum states of each sub-system, \tjj >= ^ tpia\A >i \B >a. The states \A >i are 
formed by the degrees of freedom localized in the region A while the states \B >a by those which 
are defined in the region B. The density matrix that corresponds to a pure quantum state j?/' > 

Po{A,B)^\i;><i;\ (1) 

has zero entropy. By tracing over the degrees of freedom in region A we obtain a density matrix 

PB = Tr APo{A,B) (2) 

with elements {pB)ab = {'>ljip^)ab- The statistical entropy, defined for this density matrix by the 
standard formula 

Sb = -Trps InpB (3) 

is by definition the entanglement entropy associated with the surface S. We could have traced over 
the degrees of freedom located in region B and formed the density matrix {pA)ij = (^"^V'*)ij - K 
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is clear tlia10 

Trp\ - Trp| 

for any integer k. Thus we conclude that the entropy ([3]) is the same for both density matrices pA 
and pb, 

Sa = Sb. (4) 

This property indicates that the entanglement entropy for a system in a pure quantum state is not 
an extensive quantity. In particular, it does not depend on the size of each region A ov B and thus 
is only determined by the geometry of E. 



2.2 Short-distance correlations 

On the other hand, if the entropy ([3]) is non- vanishing, this shows that in the global system there 
exist correlations across the surface E between modes which reside on different sides of the surface. 
In this review we shall consider the case when the system in question is a quantum field. The 
short-distance correlations that exist in this system have two important consequences: 

• the entanglement entropy becomes dependent on the UV cut-off e which regularizes the 
short-distance (or the large-momentum) behavior of the field system 

• to leading order in e^^ the entanglement entropy is proportional to the area of the surface S 

For a free massless scalar field the 2-point correlation function in d space-time dimensions has the 
standard form 

<c^{x),c^[y)>= , (5) 

F ~ y\ 

where fid = ^^1/2 ■ Correspondingly, the typical behavior of the entanglement entropy in d 
dimensions is 

where the exact pre-factor depends on the regularization scheme. Although the similarity between 
([5]) and ([6]) illustrates well the field-theoretical origin of the entanglement entropy, the exact relation 
between the short-distance behavior of 2-point correlation functions in the field theory and the UV 
divergence of the entropy is more subtle, as we shall discuss later in the paper. 



2.3 Thermal entropy 

Instead of a pure state one could have started with a mixed thermal state at temperature T with 
density matrix pq{A,B) = e^^ h{A,b) ^ where H{A,B) is the Hamiltonian of the global system. 
In this case the relation ^ is no more valid and the entropy depends on the size of the total 
system as well as on the size of each sub-system. By rather general arguments, in the limit of large 
volume the reduced density matrix approaches the thermal density matrix. So that in this limit 
the entanglement entropy (j3|) reproduces the thermal entropy. For further references we give here 
the expression 

Sthermal = ^r(^)C(d) T^-'V^-l (7) 

for the thermal entropy of a massless field residing inside a spatial (d— l)-volume Vd-i at temper- 
ature T. 

^For finite matrices this property indicates that the two density matrices have the same eigenvalues. 
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2.4 Entropy of a system of finite size at finite temperature 

In more general situation one starts with a system of finite size L in a mixed tliermal state 
at temperature T. Tliis system is divided by the entanghng surface S in two sub-systems of 
characteristic size I. Then the entanglement entropy is a function of several parameters (if the 
field in question is massive then mass m should be added to the parameters on which the entropy 
should depend on) 

S = Sir, L, /, e) , (8) 

where e is a UV cut-off. Clearly, the entanglement entropy in this general case is due to a combi- 
nation of different factors: the entanglement between two sub-systems and the thermal nature of 
the initial mixed state. In d dimensions even for simple geometries this function of 4 variables is 
not known explicitly. However, in two space-time dimensions, in some particular cases, the explicit 
form of this function is known. 



2.5 Entropy in (l+l)-dimensional space-time 

The state of a quantum field in two dimensions is defined on a union of intersecting intervals 
A\J B. The 2-point correlation functions behave logarithmically in the limit of coincident points. 
Correspondingly, the leading UV divergence of the entanglement entropy in two dimensions is 
logarithmic. For example, for a 2d massless conformal field theory, characterized by a central 
charge c, the entropy is |207| . [55] . [133] 

^2d = ?ln- + s(/^/ZB), (9) 
6 e 

where n is the number of intersections of intervals A and B where the sub-systems are defined, 
I A (Ib) is the length of the interval A (B). The second term in ([9]) is a UV finite term. In 
some cases the conformal symmetry in two dimensions can be used to calculate not only the 
UV divergent term in the entanglement entropy but also the UV finite term, thus obtaining the 
complete answer for the entropy, as was shown by Holzhey, Larsen and Wilczek |133| (see |161) . 
|30] for more recent developments). There are two different limiting cases when the conformal 
symmetry is helpful. In the first case one considers a pure state of the conformal field theory on 
a circle of circumference L, the subsystem is defined on a segment of size I of the circle. In the 
second situation the system is defined on an infinite line, the subsystem lives on interval of length 
/ of the line and the global system is in a thermal mixed state with temperature T. In Euclidean 
signature both geometries represent a cylinder. For a thermal state the compact direction on the 
cylinder corresponds to Euclidean time r compactified to form a circle of circumference /? = 1/T. 
In both cases the cylinder can be further conformally mapped to a plane. The invariance of the 
entanglement entropy under conformal transformation can be used to obtain 



c 



L . ,Trl 



in the case of a pure state on a circle and 



S='-ln(^shM^) (11) 



3 \T:e ^ P 

for a thermal mixed state on an infinite line. In the limit of large I the entropy (jlip approaches 
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where the first terra represents the entropy of the thermal gas ([7]) in a cavity of size / while the 
second term represent the purely entanglement contribution (note that the intersection of A and 
B contains two points in this case so that n = 2). The third term is an intermediate term due 
to the interaction of both factors, thermality and entanglement. This example clearly shows that 
for a generic thermal state the entanglement entropy is due to the combination of two factors: the 
entanglement between two subsystems and the thermal nature of the mixed state of the global 
system. 

2.6 The Euclidean path integral representation and the replica method 

A technical method very useful for the calculation of the entanglement entropy in a field theory 
is the so-called the replica trick, see ref.[33]. Here we illustrate this method for a field theory 
described by a second order Laplace type operator. One considers a quantum field ip{X) in a d- 
dimensional spacetime and chooses the Cartesian coordinates — (t, i — 1, ..,d~2), where 
r is Euclidean time, such that the surface S is defined by the condition x = and (z*, i = 1, .., d—2) 
are the coordinates on E. In the subspace (r, x) it will be convenient to choose the polar coordinate 
system r = rsin((/)) and x = rcos{(f)), where the angular coordinate (f) varies between and 2n. 
We note that if the field theory in question is relativistic then the field operator is invariant under 
the shifts — ^ 4> + w, where w is an arbitrary constant. 

One first defines the vacuum state of the quantum field in question by the path integral over 
a half of the total Euclidean spacetime defined as r < such that the quantum field satisfies the 
fixed boundary condition ijj(T = 0, x, z) = ipo{x, z) on the boundary of the half-space, 

nM^.z)]= J 2?V'e-^[^l, (13) 

where W[tlj] is the action of the field. The surface E in our case is a plane and the Cartesian 
coordinate x is orthogonal to E. The co-dimension 2 surface E defined by the conditions a; = 
and T — naturally separates the hypersurface r = in two parts: a; < and a; > 0. These are 
the two sub-regions A and B discussed in section \2A[ 

The boundary data ip(x,z) is also separated into ip-{x,z) = 'ipo{x,z), x < and Tp^ = 
'ipoix, z),x > 0. By tracing over ip- (a;, z) one defines a reduced density matrix 

p(Vi>+) - J p^-^'(^|,^-)*(v4,^-) , (14) 

where the path integral goes over fields defined on the whole Euclidean spacetime except a cut 
{t — 0,x > 0). In the path integral the field ip{X) takes the boundary value ■(/'+ above the cut 
and ip\, below the cut. The trace of the n-th power of the density matrix ([T4|) is then given by 
the Euclidean path integral over fields defined on an n-sheeted covering of the cut spacetime. In 
the polar coordinates (r, </>) the cut corresponds to values = 27rfc, A; = 1, 2, .., n. When one passes 
across the cut from one sheet to another, the fields are glued analytically. Geometrically this n-fold 
space is a flat cone C„ with angle deficit 27r(l — n) at the surface E. Thus we have 

Trp" = Z[Cn] , (15) 

where Z[Cn] is the Euclidean path integral over the n-fold cover of the Euclidean space, i.e. over 
the cone C„. Assuming that in ()15p one can analytically continue to non-integer values of n, one 
observes that 

-Trplnp = -{ada - l)lnTrp"|„=i , 
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where p — is the renormahzed matrix density. Introduce the effective action W[q\ = — \nZ{a), 

where Z{a) — Z[Ca] is the partition function of the field system in question on a Euchdean space 
with conical singularity at the surface E. In the polar coordinates (r, (p) the conical space is 
defined by making the coordinate 4> periodic with period 27ra, where (1 — a) is very small. The 
invariance under the abelian isometry (j) ^ (j) + w helps to construct without any problem the 
correlation functions with the required periodicity 27ra starting from the 27r-periodic correlation 
functions. The analytic continuation of Trp" to a different from 1 in the relativistic case is naturally 
provided by the path integral Z(a) over the conical space Ca- The entropy is then calculated by 
the replica trick 

S^iadc.-l)W{a)\^=i. (16) 

One of the advantages of this method is that we do not need to care about the normalization of 
the reduced density matrix and can deal with a matrix which is not properly normalized. 

2.7 Uniqueness of analytic continuation 

The uniqueness of the analytic continuation of Trp" to non-integer n may not seem obvious, 
especially if the field system in question is not relativistic so that there is no isometry in the polar 
angle (j) which would allow us without any trouble to glue together pieces of the Euclidean space 
to form a path integral over a conical space Ca- However, some arguments can be given that the 
analytic continuation to non-integer n is in fact unique. 

Consider a renormahzed density matrix p = t^t- The eigenvalues of p lie in the interval 
< A < 1. If this matrix were a finite matrix we could use the triangle inequality to show that 

|Trp"| < |(Trp)"| = 1 if Re{a) > 1 . 

For infinite size matrices the trace is usually infinite so that a regularization is needed. Suppose 
that e is the regularization parameter and Tr^ is the regularized trace. Then 

|Trep"| < 1 if Re{a) > 1 . (17) 

Thus Trp" is a bounded function in the complex half-plane, Re{a) > 1. Now suppose that we 
know that Ti^p^la^n = ZQ{n) for integer values of a — n, n ~ 1,2,3,... Then, in the region 
Re{a) > 1, we can represent Z(a) = Tr^p" in the form 

Z{a) ~ Zo{a) + sin(7rQ!)g(Q!) , (18) 

where the function g{a) is analytic (for Re{a) > 1). Since by condition (|17p the function Z{a) is 
bounded, we obtain that, in order to compensate for the growth of the sine in (|18p for complex 
values of a, the function g{a) should satisfy the condition 

\g{a = x + iy)\ < e^'^l^l . (19) 

By Carlson's theorem |36) an analytic function, which is bounded in the region Re(a) > 1 and 
which satisfies condition ()19|) . vanishes identically. Thus we conclude that g{a) = and there is 
only one analytic continuation to non-integer n, namely the one given by function Zo(a). 

2.8 Heat kernel and the Sommerfeld formula 

Consider for concreteness a quantum bosonic field described by a field operator T) so that the 

1/2 

partition function is Z = dot ' T>. Then the effective action defined as 

-1 /"OO 7 

W^-- -TrK{s), (20) 
2 J ^2 s 
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where parameter e is a UV cutoff, is expressed in terms of the trace of the heat kernel K(s, X, X') =< 
X\e^'''^\X' >. The latter is defined as a solution to the heat equation 

ids+V)Kis,X,X') = 0, 

K{s^O,X,X') = SiX,X'). ^ ' 

In order to calculate the effective action Vl^(a) we use the heat kernel method. In the context of 
manifolds with conical singularities this method was developed in great detail in [69], [103) . In the 
Lorentz invariant case the invariance under the abelian symmetry ^ — >■ + w plays an important 
role. The heat kernel K{s, (j), (j)') (where we omit the coordinates other than the angle (f)) on regular 
flat space then depends on the difference (0 — </>'). This function is 27r periodic with respect to 
{(j) — (j)'). The heat kernel Ka{s, (j), (j)') on a space with a conical singularity is supposed to be 27ra 
periodic. It is constructed from the 27r periodic quantity by applying the Sommerfeld formula |206j 

Z f ID 

Ka{s, (t), (/)') = K{s,(f>- (f)') + / cot —K{s, - </>' + w)dw . (22) 

Ana Jp 2a 

That this quantity still satisfies the heat kernel equation is a consequence of the invariance under 
the abelian isometry — ^ (j) + w. The contour F consists of two vertical lines, going from (— tt + ioo) 
to (— TT — ioo) and from (tt — ioo) to (tt — |-«oo) and intersecting the real axis between the poles of 
the cot — 27ra, and 0, +27ra, respectively. For a = 1 the integrand in ([22]) is a 27r-periodic 
function and the contributions of these two vertical lines cancel each other. Thus, for a small angle 
deficit the contribution of the integral in (|22p is proportional to {1 — a). 

2.9 An explicit calculation 

Consider an infinite (d— 2)-plane in d-dimensional space-time. The calculation of the entanglement 
entropy for this plane can be done explicitly by means of the heat kernel method. In fiat space-time, 
if the operator V is the Laplace operator, 

V = -V^ 

one can use the Fourier transform in order to solve the heat equation. In d spacetime dimensions 
one has 

K{s,X,X') = I d'^pe^f-(^"-^'") e-'^^P"'> . (23) 

(Stt)'* j 

Putting z' = z'*, i = 2 and choosing in the polar coordinate system {r^cj)), that <j) = (j)' + w 

we have that p^{X—X'Y — 2pr sin ^ cos 0, where = P^Pfi and 9 is the angle between the d- 
vectors and {Xf^—X'^). The radial momentum p and angle 9, together with the other (rf— 2) 
angles form a spherical coordinate system in the space of momenta p^^. Thus one has for the 
integration measure / d'^p = il.d-2 /g°° dpp'^^'^ d9 sin''^^ 9 , where ftd-2 = r({d-i)/2) ^^^^ 
of a unit radius sphere in d—1 dimensions. Performing the integration in (j23p in this coordinate 
system we find 

^^^^JZ^ I dpp^J._^{2rps.n-)e - . (24) 



For the trace one finds 



TrK{s,w) = / .T^ A(S) , (25) 
(47rs)2 sm^ f 
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where ^(S) — J d'^ "^z is the area of the surface S. One uses the integral dxx^ ^ J^ix) = \jp^ 
for the derivation of (psj) . The integral over the contour F in the Sommerfeld formula (|22p is 
calculated via residues ([59], |103J ) 

C2 a =— / cot— = l-a^ . 26 
Collecting everything together one finds that in flat Minkowski space-time 

^ (4^)^/2 ("^ + 2™C2(a)sA(I])) , (27) 

where V — J drd'^^^x is the volume of space-time and A{Y^) = J d'^^'^x is the area of the surface 
S. Substituting (|27p into equation ([20|) we obtain that the effective action contains two terms. 
The one proportional to the volume V reproduces the vacuum energy in the effective action. The 
second term proportional to the area ACS) is responsible for the entropy. Applying the formula 
(jl6p we obtain the entanglement entropy 

6(d-2)(4^)(rf-2)/2g<i-2 ^ J 

of an infinite plane E in d space-time dimensions. Since any surface, locally, looks like a plane 
and a curved spacetime, locally, is approximated by Minkowski space, this result gives the leading 
contribution to the entanglement entropy of any surface S in flat or curved space-time. 

2.10 Entropy of massive fields 

The heat kernel of a massive field described by the wave operator P = — + is expressed in 
terms of the heat kernel of massless field, 

K(m^O) (X, x' , S) = K(^rn=Q) (x, x' , s) • 6^" \ 

Thus one finds 

Trif^™^°) (s) = TrX^™=°) (s) • e'""'' , (29) 

where the trace of the heat kernel for vanishing mass is given by (|27|). Therefore the entanglement 
entropy of a massive field is 

_ Ajj:) ds 

" 12(47r)('^-2)/2 ^dT^^ • i^U) 
In particular, if d = 4, one finds that 
A(T,) 1 

S^^o = —}-^{^ + 2m^lne + m^lnm''+m^{-f-l) + 0{em)). (31) 
12(47r) 

The logarithmic term in the entropy that is due to the mass of the field appears in any even 
dimension d. The presence of a UV finite term proportional to the (d — 2)-th power of mass is the 
other general feature of ([50]) . (|5T|) . 
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2.11 An expression in terms of the determinant of the Laplacian on the 
surface 



Even though the entanglement entropy is determined by the geometry of the surface S, in general 
this can be not only its intrinsic geometry but also how the surface is embedded in the larger 
space-time. The embedding is determined by the extrinsic curvature. The curvature of the larger 
spacetime enters through the Gauss-Cadazzi relations. But in some particularly simple cases the 
entropy can be given a purely intrinsic interpretation. To see this for the case when E is a plane we 
note that the entropy (pS)) or (I30p originates from the surface term in the trace of the heat kernel 
((77)) (or ([5T|)). To leading order in (1 — a), the surface term in the case of a massive scalar field is 

(l_a).i.Tri^s(s), 



where 

Tri^s(s) 



A(S) 



(47rs) 2 

can be interpreted as the trace of the heat kernel of operator — A(E) + rn^, where A(E) is the 
intrinsic Laplace operator defined on (d— 2)-plane E. The determinant of the operator — A(E) -|-to^ 
is determined by 

P°° ds 

lndet(-A(I]) + m^) = - / —TrKY:{s) . 



s 

Thus we obtain an interesting expression for the entanglement entropy 

S' = -^lndet(-A(I])+m2) (32) 

in terms of geometric objects defined intrinsically on the surface E. A similar expression in the 
case of an ultra-extreme black hole was obtained in |172j and for a generic black hole with horizon 
approximated by a plane was obtained in [97] . 



2.12 Entropy in theories with a modified propagator 

In certain physically interesting situations the propagator of a quantum field is different from the 
standard l/{p'^ + 'm?) and is described by some function as 1/F(p'^). The quantum field in question 
then satisfies a modified Lorentz invariant field equation 

= F(V2)V' = 0. (33) 

Theories of this type naturally arise in models with extra dimensions. The deviations from the 
standard form of propagator may be both in the UV regime (large values of p) or in the IR regime 
(small values of p). If the function F{p'^) for large values of p grows faster than p^ this theory is 
characterized by improved UV behavior. 

The calculation of the entanglement entropy performed in section 12.91 can be generalized to 
include theories with operator p3p . This example is instructive since, in particular, it illuminates 
the exact relation between the structure of 2-point function (the Green's function in the case of 
free fields) and the entanglement entropy |182J . 

In d spacetime dimensions one has 

K{s,X,Y) = d-^pe'^^^^"-^^) e'^^(^=) . (34) 
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Note that we consider Euclidean theory so that > 0. The Green's function 

G(x,y)=<7A(x),v(r)> (35) 

is a solution to the field equation with a delta-like source 

VGiX,Y) ^ S{X,Y) (36) 
and can be expressed in terms of the heat kernel as follows 



/■oo 

G{X,Y)= / dsK{s,X,Y) . (37) 
Jo 

Obviously, the Green's function can be represented in terms of the Fourier transform in a manner 
similar to p4|) . 

G{X, Y) = ^ j d'p e'f-(^^-^") G(p2) , G(p2) = i/F{p') . (38) 

The calculation of the trace of the heat kernel for operator (|33l) on a space with a conical singularity 
goes along the same lines as in section [^751 This was performed in |183] and the result is 

Tri^„(s) - ^^-L_(aFFrf(s) + 2™G2(a)^(S)Prf-2(s)) , (39) 

where the functions Pn{s) are defined as 

Pn{^) = Y(^l dpp"-'e-^^^^'K (40) 

The entanglement entropy takes the form (we remind the reader that for simplicity we take the 
surface S to be a (d — 2)-dimensional plane) |183j 

s- .. r-p^~2is), (41) 



12 • (47r)('^-2)/2 
It is important to note that |183| . [182] 

(i) the area law in the entanglement entropy is universal and is valid for any function F{p^)] 

(ii) the entanglement entropy is UV divergent independently of the function F{p'^), with the degree 
of divergence depending on the particular function F{p'^); 

(iii) in the coincidence limit, X = Y, the Green's function ([38|) 

(42) 



9 1 1"^ 

G(^'^)-^77Tt/ dpp'^'Gip^) 
1 ( 9 ) 47r 2 Jo 



Til) (47r) 

may take a finite value if G{p^) = 1/F{p'^) is decaying faster than l/p"^. However, even for this 
function F(p^) the entanglement entropy is UV divergent. 

As an example, consider a function which grows for large values of p as 
2-point correlation function in this theory behaves as 

<HX),HY) >^ ,^_y,,,2. (43) 
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and for k > d/2 it is regular in the coincidence limit. On the other hand, the entanglement entropy 
scales as 



s^m («) 

e 

and remains divergent for any positive value of k. Comparison of ()43|) and (|44p shows that only for 
fc = 1 (the standard form of the wave operator and the propagator) the short-distance behavior of 
the 2-point function is similar to the UV divergence of the entanglement entropy. 

2.13 Entanglement entropy in non-Lorentz invariant theories 

Non-Lorentz invariant theories are characterized by a modified dispersion relation, uj'^ + F{p'^) = 0, 
between the energy uj and the 3-momentum p. These theories can be described by a wave operator 
of the following type 

V = -df + Fi-A,) , (45) 

where A^^ = ^^e spatial Laplace operator. Clearly, the symmetry with respect to the 

Lorentz boosts is broken in operator (|45p if F{q) ^ q. 

As in the Lorentz invariant case to compute the entanglement entropy associated with a surface 
E we choose {d — 1) spatial coordinates {x*, « = 1, .., d — 1} = {x, z"" , a — 1, ..,d — 2}, where x 
is the coordinate orthogonal to the surface S and z'^ are the coordinates on the surface S. Then, 
after going to Euclidean time r — it, we switch to the polar coordinates, t — r sin((/)), x — r cos((/)). 
In the Lorentz invariant case the conical space which is needed for calculation of the entanglement 
entropy is obtained by making the angular coordinate (p periodic with period 2'Ka by applying 
the Sommerfeld formula (|22|) to the heat kernel. If Lorentz invariance is broken, as it is for the 
operator (|45p there are certain difficulties in applying the method of the conical singularity when 
one computes the entanglement entropy. The difficulties come from the fact that the wave operator 
2?, if written in terms of the polar coordinates r and 0, becomes an explicit function of the angular 
coordinate cf). As a result of this, the operator T) is not invariant under shifts of 4> to arbitrary 
(j) + w. Only shifts with w — 2?™, where n is an integer are allowed. Thus, in this case one cannot 
apply the Sommerfeld formula since it explicitly uses the symmetry of the differential operator 
under shifts of angle On the other hand, a conical space with angle deficit 27r(l — n) is exactly 
what we need to compute Trp" for the reduced density matrix. In reHl83 , by using some scaling 
arguments it was shown that the trace of the heat kernel K{s) — e^''^ on a conical space with 
2?™ periodicity, is 

TrX„(s) = nTrif„=i(s) + .^^-L_27rnC2(n) A(E) Pd-2{s) , (46) 

where nTri<'„=i(s) is the bulk contribution. By the arguments presented in Section [2771 there is a 
unique analytic extension of this formula to non-integer n. A simple comparison with the surface 
term in the heat kernel of the Lorentz invariant operator, which was obtained in the previous 
section, shows that the surface terms of the two kernels are identical. We thus conclude that the 
entanglement entropy is given by the same formula 

where Pn{s) is defined in (HO]) , as in the Lorentz invariant case (|4ip . A similar property of the 
entanglement entropy was observed for a non-relativistic theory described by the Schrodinger 
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operator j204J (see also [59] for a holographic derivation). For polynomial operators, F{q) ~ q'', 
some scaling arguments can be used [204] to get the form of the entropy that follows from (|T7)) . 

In the rest of the review we shall mostly focus on the study of Lorcntz invariant theories with 
field operator quadratic in derivatives, of the Laplace type, V = — (V^ + X). 

2.14 Arbitrary surface in curved space-time: general structure of UV 
divergences 

The definition of the entanglement entropy and the procedure for its calculation generalize to 
curved spacetime. The surface E can then be any smooth closed co-dimension two surfaceH which 
divides the space in two sub-regions. In the next section we will consider in detail the case where 
this surface is a black hole horizon. Before proceeding to the black hole case we would like to 
specify the general structure of UV divergent terms in the entanglement entropy. In d-dimensional 
curved spacetime entanglement entropy is presented in the form of a Laurent series with respect 
to the UV cutoff e (for d = 4 see [205] ) 

„ Sd-2 , Sd-4 , , Sd-2-2n , , i , / \ fAn\ 

^ - ^ + ^ + - + gd-2-2» + - + go Ine + s{g) , (48) 

where Sd-2 is proportional to the area of the surface E. All other terms in the expansion (|48p can 
be presented as integrals over E of local quantities constructed in terms of the Riemann curvature 
of the spacetime and the extrinsic curvature of the surface E. The intrinsic curvature of the surface 
E of course can be expressed in terms of TZ and k using the Gauss-Codazzi equations. Since nothing 
should depend on the direction of vectors normal to E, the integrands in expansion pS]) should 
be even powers of extrinsic curvature. The general form of the Sd-2-2n term can be symbolically 
presented in the form 

l+p—n ' 



E 

where TZ stands for components of the Riemann tensor and their projections onto the sub-space 
orthogonal to E and k labels the components of the extrinsic curvature. Thus, since the integrands 
are even in derivatives, only terms fr''-2n-2^ n = 0, 1, 2, .. appear in (|48p . If d is even then there 
also may appear a logarithmic term sq. The term s{g) in (|48p is a UV finite term, which also may 
depend on the geometry of the surface E, as well as on the geometry of the space-time itself. 



3 Entanglement entropy of non-degenerate Killing horizons 

3.1 The geometric setting of black hole spacetimes 

The notion of entanglement entropy is naturally applicable to a black hole. In fact, probably the 
only way to separate a system in two sub-systems is to place one of them inside a black hole 
horizon. The important feature that, in fact, defines the black hole is the existence of a horizon. 
Many useful definitions of a horizon are known. In the present paper we shall consider only the 
case of the so-called eternal black holes for which different definitions of the horizon coincide. 
The corresponding spacetime then admits a maximal analytic extension which we shall use in 
our construction. The simplest example is the Schwarzschild black hole, the maximal extension 
of which is demonstrated on the well known Penrose diagram. The horizon of the Schwarzschild 

■^If the boundary of S is not empty there could be extra terms in the entropy proportional to the "area" of the 
boundary 9S as was shown in |108| . We do not consider this case here. 
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black hole is an example of a so-called Killing horizon. The spacetime in this case possesses a 
global Killing vector, ^4 = dt, which generates time translations. The Killing horizon is defined as 
a null hypersurface on which the Killing vector is null, = 0. The null surface in the maximal 
extension of an eternal black consists of two parts: the future horizon and the past horizon. The two 
intersect on a co-dimension two compact surface S called the bifurcation surface. In the maximally 
extended spacetime a hypersurface Ht of constant time i is a Cauchy surface. The bifurcation 
surface E naturally splits the Cauchy surface in two parts, and "H^, respectively inside and 
outside the black hole. For asymptotically flat spacetime, a such as the Schwarzschild metric, the 
hypersurface fit has the topology of a wormhole. (In the case of the Schwarzschild metric it is 
called the Einstein- Rosen bridge.) The surface E is the surface of minimal area in Ht- In fact the 
bifurcation surface S is a minimal surface not only in the {d — l)-dimensional Euclidean space Ht 
but also in the d-dimensional spacetime. As a consequence, as we show below, the components of 
the extrinsic curvature defined for two vectors normal to S, vanish on S. 

The space-time in question admits a Euclidean version by analytic continuation t ir. It 
is a feature of regular metrics with a Killing horizon that the direction of Euclidean time t is 
compact with period 2t:/3h which is determined by the condition of regularity, i.e. the absence of 
a conical singularity. In the vicinity of the bifurcation surface E, the spacetime then is a product 
of a compact surface E and a two dimensional disk, the time coordinate r playing the role of the 
angular coordinate on the disk. The latter can be made more precise by introducing a new angular 
variable (f> — PJj^t which varies from to 27r. We consider the static space-time with Euclidean 
metric of the general type 

ds^ = plg{p)d(t)'^ + + ^^^(^p^ e)d9'de^ . (50) 

The radial coordinate p is such that the surface E is defined by the condition p — Q. Near this 
point the functions g{p) and Jijip, d) can be expended as 

9{p) ^ i + 0{p') , j,,(p,9)=j§\e)+0{p'), (51) 
Ph 

where is the metric on the bifurcation surface E equipped with coordinates {9^ , i = 1, .., d— 

2}. This metric describes what is called a non-degenerate horizon. The Hawking temperature of 
the horizon is finite in this case and equal to Th = l/(27r /?//). 

It is important to note that the metric (|50p does not have to satisfy any field equations. The 
entanglement entropy can be defined for any metric which possesses a Killing type horizon. In 
this sense the entanglement entropy is an off-shell quantity. It is useful to keep this in mind 
when one compares the entanglement entropy with some other approaches in which an entropy is 
assigned to a black hole horizon. Even though the metric (|50l) with (|51l) does not have to satisfy 
the Einstein equations we shall still call the complete space described by the Euclidean metric dSH]) 
the Euclidean black hole instanton and will denote it by E. 

3.2 Extrinsic curvature of horizon, horizon as a minimal surface 

The horizon surface E defined by the condition p = in the metric (|50p is a co-dimension 2 surface. 
It has two normal vectors: a spacelike vector with the only non-vanishing component = 1 
and a timelike vector n'^ with the non- vanishing component n| = 1/ p. With respect to each normal 
vector one defines an extrinsic curvature, kfj = —Ji''Jj^'^inp, a = 1,2. The extrinsic curvature 
kfj identically vanishes. It is a consequence of the fact that is a Killing vector which generates 
the time translations. Indeed, the extrinsic curvature can be also written as a Lie derivative, 
fcpi/ = — ^>C„g^jy, so that it vanishes if rt is a Killing vector. The extrinsic curvature associated to 
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the vector n , 

fci- =-^7.S/5p7fc«, (52) 

is vanishing when restricted to the surface defined by the condition p = 0. It is due to the fact that 
the term hnear in p is absent in the p-expansion for ^ij{p, 0) in the metric (|50|) . This is required by 
the regularity of the metric (|50p : in the presence of such a term the Ricci scalar would be singular 
at the horizon, R ^ ^1 P- 

The vanishing of the extrinsic curvature of the horizon indicates that the horizon is necessarily 
a minimal surface. It has the minimal area considered as a surface in d-dimensional spacetime. On 
the other hand, in the Lorentzian signature, the horizon S has the minimal area if considered on 
the hypersurface of constant time t, Ht, the latter thus has the topology of a wormhole. 

3.3 The wave function of a black hole 

Although the entanglement entropy can be defined for any co-dimension two surface in spacetime 
when the surface is a horizon particular care is required. In order to apply the general prescription 
outlined in section 12. 1[ we first of all need to specify the corresponding wave function. Here we 
will follow the prescription proposed by Barvinsky, Frolov and Zelnikov [16j . This prescription is 
a natural generalization of the one in flat spacetime discussed in section [Z6] On the other hand, 
it is similar to the "no-boundary" wave function of the Universe introduced in [127' . We define 
the wave function of a black hole by the Euclidean path integral over field configurations on the 
half-period Euclidean instanton defined by the metric (|50p with angular coordinate (j) changing 
in the interval from to it. This half-period instanton has Cauchy surface H (on which we can 
choose coordinates x = {p, 9)) as a boundary where we specify the boundary conditions in the path 
integral, 

*[7/'_(x),V+(a;)] = J 2?Ve-'^['^l, (53) 

i:{X)\^^o ^ Mx) 
i;{X)\^^^ ^ yj^ix) 

where Wltp] = \ j tpVtp is the action of the quantum field V'- The functions tp-{x) and tlj+{x) are 
the boundary values defined on the part of the hypersurface "H which is respectively inside (T-l-) 
and outside (71+) the horizon S. As was shown in [16] the wave function ([53)1 corresponds to the 
Hartle-Hawking vacuum state |126j . 

3.4 Reduced density matrix and entropy 

The density matrix p('!/'+,V'+) defined by tracing over ip^-modes is given by the Euclidean path 
integral over field configurations on the complete instanton {0 < (p < 2tt) with a cut along the axis 
</) = where the field V'(A^) in the path integral takes the values ■0+(2^) and ip\{x) below and above 
the cut respectively. The trace Trp is obtained by equating the fields across the cut and doing the 
unrestricted Euclidean path integral on the complete Euclidean instanton E. Analogously, Trp" is 
given by the path integral over field configurations defined on the n-fold cover £"„ of the complete 
instanton. This space is described by the metric (|50p where angular coordinate cf) is periodic with 
period 2?™. It has a conical singularity on the surface E so that in a small vicinity of S the total 
space En is a direct product of E and a two-dimensional cone C„ with angle deficit 5 = 27r(l — n). 
Due to the abelian isometry generated by the Killing vector this construction can be analytically 
continued to arbitrary (non-integer) n — > a. So that one can define a partition function 

Z{a) = Trp" (54) 
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by the path integral over field configurations over Ea, the a- fold cover of the instanton E. For 
a bosonic field described by the field operator V one has that Z{a) — det^^^^V. Defining the 
effective action as VF(a) = — InZ(a), the entanglement entropy is still given by (|16l) . i.e. by 
differentiating the effective action with respect to the angle deficit. Clearly, only the term linear 
in (1 — a) contributes to the entropy. The problem thus reduces to the calculation of this term in 
the effective action. 

3.5 The role of the rotational symmetry 

We emphasize that the presence of the so called rotational symmetry with respect to the Killing 
vector 9^, which generates rotations in the 2-plane orthogonal to the entangling surface S, plays 
an important role in our construction. Indeed, without such a symmetry it would be impossible 
to interpret Trp" for an arbitrary a as a partition function in some gravitational background. In 
general, two points are important for this interpretation: 

i) that the spacetime possesses, at least locally near the entangling surface, a rotational symmetry 
so that, after the identification — (/) + 27ra we get a well-defined spacetime E^ with no more than 
just a conical singularity; this holds automatically if the surface in question is a Killing horizon; 

ii) and that the field operator is invariant under the "rotations", — >■ + w; this is automatic if 
the field operator is a covariant operator. 

In particular, the point ii) allows us to use the Sommerfeld formula (more precisely its gener- 
alization to a curved spacetime) in order to define the Green's function or the heat kernel on the 
space Ea- As is shown in 183 (see also discussion in section [2. 131) in the case of the non-Lorentz 
invariant field operators in fiat Minkowski spacetime the lack of the symmetry ii) makes the whole 
"conical space" approach rather obscure. On the other hand, in the absence of the rotational sym- 
metry i) there may appear terms in the entropy that are "missing" in the naively applied conical 
space approach: the extrinsic curvature contributions [203] or even some curvature terms [134] . 

In what follows we consider the entanglement entropy of the Killing horizons and deal with the 
covariant operators so that we do not have to worry about i) or ii). 

3.6 Thermality of the reduced density matrix of a Killing horizon 

The quantum state defined by equation ([53]) is the Hartle- Hawking vacuum [126) . The Green's 
function in this state is defined by analytic continuation from the Euclidean Green's function. The 
periodicity t — t + iPn is thus inherent in this state. This periodicity indicates that the correlation 
functions computed in this state are in fact thermal correlation functions when continued to the 
Lorentzian section. This fact generalizes to an arbitrary interacting quantum field as shown in 
|121) . On the other hand, being globally defined, the Hartle-Hawking state is a pure state which 
involves correlations between modes localized on different sides of the horizon. This state however 
is described by a thermal density matrix if reduced to modes defined on one side of the horizon 
as was shown by Israel |138j . That the reduced density matrix obtained by tracing over modes 
inside the horizon is thermal can be formally seen by using angular quantization. Introducing 
the Euclidean Hamiltonian He which is the generator of rotations with respect to the angular 
coordinate defined above, one finds that p(^+, ')/'+) =< IV'+ >j i-S- the density matrix 

is thermal with respect to the Hamiltonian He with inverse temperature 27r. This formal proof 
in Minkowski space was outlined in |152j . The appropriate Euclidean Hamiltonian is then the 
Rindler Hamiltonian which generates the Lorentz boosts in the direction orthogonal to the surface 
S. In [141] the proof was generalized to the case of generic static spacetimes with bifurcate Killing 
horizons admitting a regular Euclidean section. 
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3.7 Useful mathematical tools 



3.7.1 Curvature of space with a conical singularity 

Consider a space Ea which is an a-fold covering of a smooth manifold E along the Killing vector 
dip generating an Abelian isometry. Let surface S be a stationary point of this isometry so that 
near E the space Ea looks like a direct product, E x C^, of the surface E and a two-dimensional 
cone Ca with angle deficit S — 2tt{1 — a). Outside the singular surface E the space Ea has the 
same geometry as a smooth manifold E. In particular, their curvature tensors coincide. However, 
the conical singularity at the surface E produces a singular (delta-function like) contribution to 
the curvatures. This was first demonstrated by Sokolov and Starobinsky [194] in two-dimensions 
by using topological arguments. These arguments were generalized to higher dimensions in [7]. 
One way to extract the singular contribution is to use some regularization procedure, replacing 
the singular space Ea by a sequence of regular manifolds Ea- This procedure was developed by 
Fursaev and Solodukhin in |112| . In the limit Ea Ea one obtains the following results |112J : 

^""a/J - ^''1/3 + 27r(l - a) ({n'^na){n''n^) - {n^np){n''na)) , 
Rt^^ = Rt^^ + 27r(l - a)(n^n^)(5E , 

i? = i? + 47r(l - a)(SE , (55) 

where Ss is the delta-function, fds ~ /j. /; n'^ — n'^d^ , k — 1, 2 are two orthonormal vectors 

orthogonal to the surface E, (n^nj^) = "^fi'^^ ^^'^ quantities R^'^ap^ ^^'^ ^ 

computed in the regular points Eaj^ by the standard method. 

These formulas can be used to define the integral expressiontH [112] 

/ R = a [ R + 4:TT{l-a) ( 1 , (56) 

R^ = a I R^ + 8TT{l~a) [ R + 0{{l-af ) , (57) 
Je jt, 



R^'^R^.u = a / + 4^(1 - a) / Ru + 0((1 - aY) , (58) 

J E JT. 

[ R'^'^PR^.^xp = a [ Rf^'^PR^^xp + 8^(1 - a) [ + 0((1 - af) , (59) 
JEa Je Jt. 

where Ru — -R^j^nf nj' and Rijij = Rfj,,^\pn^n^njnP. We use a shorthand notation for the surface 
integral = ^d''~^e. 

The terms proportional to a in (|56l) - (|59l) are defined on the regular space E. The terms 
0((1 — a)^) in (|57|) - (|59|) are something like a square of the (5-function. They are not well-defined 
and depend on the way the singular limit Ep S/j is taken. Those terms however are not important 
in the calculation of the entropy since they are of higher order in (1 — a) . There are however certain 
invariants, polynomial in the Riemann tensor, in which the terms 0((1 — a)^) do not appear at 
all. These invariants are thus well defined on the manifolds with conical singularity. Below we 
consider two examples of such invariants |112) . 



^It should be noted that formulas Il55| l and l|56ll . Il57| l. I|58ll . I|59|l are valid even if subleading terms (as in l|51ll ) 
in the expansion of the metric near singular surface S are functions of 6 11121 . Such more general metrics describe 
what might be called a "local Killing horizon" . 
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Topological Euler number. The topological Eulcr number of a 2p-diniensional smooth mani- 
fold £ is given by the integral 

J E 

C ~ n f ^VW-l.-.V^-p-W^p Tjllllll TDfJ.2p-lfJ.2p r 

02p — '^p'^fJ.lfl2---fJ.2p-lH2p'^ ^ Uiiy2---^ l^2p-11^2p ' 

Suppose that Ea has several singular surfaces (of dimension 2{p — 1)) S.^, each with conical deficit 
27r(l — ai). The Euler characteristic of this manifold is jll2l 

xiEc] = / £2p + V(l - . (61) 

Je^/s 

A special case is when Ea possesses a continuous Abelian isometry. The singular surfaces Si are 
the fixed point sets of this isometry so that all surfaces have the same angle deficit ai = a. The 
Euler number in this case is [112J 

x[Ea] = axlEa^i] + {l-a)Y^ x[^^] ■ (62) 

i 

An interesting consequence of this formula is worth mentioning. Since the introduction of a conical 
singularity can be considered as the limit of certain smooth deformation, under which the topolog- 
ical number does not change, one has — x[Ea=i]- Then one obtains an interesting formula 
reducing the number x of a manifold E to that of the fixed points set of its abelian isometry |112| 

X[i?a=i] =^Xp.]. (63) 

i 

A simple check shows that (|63|) gives the correct result for the Euler number of the sphere S^- 
Indeed, the fixed points of 2-sphere are its "north" and "south" poles. Each of these points 
has X = 1 and one gets from (|63p : xl^"^] = 1 + 1 = 2. On the other hand, the singular surface of 
{d > 3) is S"^-2 and from §^ the known identity x[5''^] = x[5"'"^] follows. Note that equation 
l|63p is valid for spaces with continuous abelian isometry and it may be violated for an orbifold 
with conical singularities. 



2^PTTPp 



^ (60) 



Lovelock gravitational action. The general Lovelock gravitational action is introduced on a 
d-dimensional Riemannian manifold as the following polynomial 166 

p=i p=i 

where ' j is the totally antisymmetrized product of the Kronecker symbols and fc^ is (d— 2)/2 
(or [d — l)/2) for even (odd) dimension d. If the dimension of spacetime is 2p, the action 
reduces to the Euler number (|60l) and is thus topological. In other dimensions the action (|64|) is 
not topological, although it has some nice properties which make it interesting. In particular, the 
field equations which follow from (j64l) are quadratic in derivatives even though the action itself is 
polynomial in curvature. 

*Note that in ref. 11121 there is a typo in eq.(3.9) defining Cp. This does not affect the conciusions of ref. I112| 
since they are based on the relation Cp_i = SirpCp rather than on the explicit form of Cp. 
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On a conical manifold Aia, the Lovelock action is the sum of volume and surface parts |112] 
Wl[Mo.] = Wl[Mo./^ + 27r(l -a)Y, Ap+iW^p[S] , (65) 

p=0 

where the first term is the action computed at the regular points. As in the case of the topological 
Euler number, all terms quadratic in (1 — a) mutually cancel in ()65p . The surface term in ()65|) 
takes the form of the Lovelock action on the singular surface S. It should be stressed that integrals 
Wp [S] are defined completely in terms of the intrinsic Riemann curvature R''-' of E 

and Wq = /j,. Ea. (|55|) allows one to compute the entropy in the Lovelock gravity by applying 
the replica formula. In ref. [145] this entropy was derived in the Hamiltonian approach, whereas 
arguments based on the dimensional continuation of the Euler characteristics have been used for 
its derivation in [7]. 



3.7.2 The heat kernel expansion on a space with a conical singularity 

The useful tool to compute the effective action on a space with a conical singularity is the heat 
kernel method already discussed in section 12.81 In section 12.91 we have shown how, in flat space, 
using the Sommerfeld formula (|22l) . to compute the contribution to the heat kernel due to the 
singular surface S. This calculation can be generalized to an arbitrary curved space Ea that 
possesses, at least locally, an Abelian isometry with a flxed point. To be more specific we consider 
a scalar field operator V — — (V^ + X), where X is some scalar function. Then the trace of the 
heat kernel K = e"''^ has the following small s expansion 

Trif£„(s) = -^5]a„s", (67) 

where the coefficients in the expansion decompose into bulk (regular) and surface (singular) parts 

a„ = a:^'s + . (68) 

The regular coefficients are the same as for a smooth space. The first few coefficients are 

The coefficients due to the singular surface E (the stationary point of the isometry) are 

3 a 

^_7r (l-a)(l + a) f I n (1 ~ a)(l + a)(l + g^) f 



The form of the regular coefficients (15^ in the heat kernel expansion has been well studied in 
physics and mathematics literature (for a review see |218| ). The surface coefficient of in (170)) was 
calculated by the mathematicians McKean and Singer |174) (see also |42^). In physics literature 
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this term has appeared in the work of Dowker [69]. (In the context of cosmic strings one has 
focused more on the Green's function rather on the heat kernel [3], llOOj .) The coefficient has 
been first obtained by Fursaev [103] although in some special cases it was known before in works 
of Donnelly [64], [65]. 

It should be noted, that due to the fact that the surface S is a fixed point of the Abelian 
isometry, all components of the extrinsic curvature of the surface S vanish. This explains why the 
extrinsic curvature does not appear in the surface terms ()70|) in the heat kernel expansion. 



3.8 General formula for entropy in the replica method, relation to the 
Wald entropy 

As a consequence of the expressions ([55]) for the curvature of space with a conical singularity 
that were presented in section [3.7.11 one obtains a general expression for the entropy. Consider a 
Euclidean general covariant action 

W[g^,,^A]^~J d''x^£{g^,,R''l,,W,R''l„...,^A), (71) 

which describes the gravitational field coupled to some matter fields ipA- In the replica trick we 
first introduce a conical singularity at the horizon surface S with a small angle deficit 5 = 27r(l — a) 
so that the Riemann curvature obtains a delta-like surface contribution (|55|1 and the gravitational 
action (|7ip becomes a function of a. Then applying the replica formula 

S ={ada,-l)W{a)\^^i 

we get 

S = 2n [ Q^p^, ((7^^n")(n"^n^) - (n'^n'5)(n''n")) (72) 



for the entropy associated to S, where tensor Qap^iu is defined as a variation of the action ()7ip 
with respect to the Riemann tensor. 



(73) 



If the action (fTT]) is local and it does not contain covariant derivatives of the Riemann tensor then 



the tensor is a partial derivative of the Lagrangian, 



Now, as was observed by Myers and Sinha |180] (see also [4]), one can re-express 

2 

^ («)(n>,^) - «r.f )(n>p = , (75) 

where e"^ — riin^ — n'^n^ is the two-dimensional volume form in the space transverse to the 
horizon surface S. Then for a local action (|71l) polynomial in the Riemann curvature the entropy 
(|72]) takes the form 

r FjC 

T^^M.^"^ , (76) 
JT, on 7,1/ 
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which is exactly the Wald entropy j220) . [144) . It should be noted that Wald's Noether charge 
method is an on-shell method so that the metric in the expression for the Wald entropy is supposed 
to satisfy the field equations. On the other hand, the conical singularity method is an off-shell 
method valid for any metric that describes a black hole horizon. The relation between the on-shell 
and the off-shell descriptions will be discussed in section 



3.9 UV divergences of entanglement entropy for a scalar field 

For a bosonic field described by a field operator 2? the partition function is Z{a) — det'^^^V. 
The corresponding effective action VK(a) = — lnZ{a) on a space with a conical singularity, E^, is 
expressed in terms of the heat kernel Ke^ (s) in a standard way 

1 A 

W{a) = -- -Trif^Js), (77) 

2 7,2 s 

The entanglement entropy is computed using the replica trick as 

S=iad^-l)Wia)\^^i. (78) 

Using the small s expansion one can, in principle, compute all UV divergent terms in the entropy. 
The surface terms are however known only for the first few terms in the expansion (|67p. This 
allows us to derive an explicit form for the UV divergent terms in the entropy. 



In two dimensions the horizon is just a point and the entanglement entropy diverges logarith- 
mically [33], [152], [70], [85], [l97] 



Sd=2 = 7 In - 
6 e 



(79) 



In three dimensions the horizon is a circle and the entropy 



^.^3 = (80) 
12-y/7re 



is linearly divergent. 



The leading UV divergence in d dimensions can be computed directly by using the form 
of the coefficient af l|70p in the heat kernel expansion [33] 

^'=6(d-2)'(4.)^^^- ^''^ 

It is identical to expression ()28|) for the entanglement entropy in flat Minkowski spacetime. This has 
a simple explanation. To leading order the spacetime near the black hole horizon is approximated 
by the fiat Rindler metric. The leading UV divergent term in the entropy is thus the entanglement 
entropy of the Rindler horizon. The curvature corrections then show up in the subleading UV 
divergent terms and in the UV finite terms. 
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The four-dimensional case is the most interesting since in this dimension there appears a 
logarithmic subleading term in the entropy. For a scalar field described by a field operator — (V^ + 
X) the UV divergent terms in the entanglement entropy of a generic 4-dimensional black hole read 

m 

We note that for a massive scalar field X — —m^. 

Of special interest is the case of the 4d conformal scalar field. In this case X — —\R and the 
entropy (|82p takes the form 

The logarithmic term in (|83|) is invariant under the simultaneous conformal transformations of 
bulk metric g^j^^, — > e^" g^i, and the metric on the surface E, — e^'^^ij. This is a general feature 
of the logarithmic term in the entanglement entropy of a conformally invariant field. 

Let us consider some particular examples. 



3.9.1 The Reissner-Nordstrom black hole. 

A black hole of particular interest is the charged black hole described by the Reissner-Nordstrom 
metric, 



ds|j^ = g{r)dT^ + g \r)dr^ + r\d6^ + sin' 
g{r) = l-^^-I±)^^. (84) 

This metric has a vanishing Ricci scalar, R — 0. It has inner and out horizons, r_ and r+ 
respectively, defined by 



r± 



. ± ^m? - g2 , (85) 



where m is the mass of the black hole and q is the electric charge of the black hole. The two 
vectors normal to the horizon are characterized by the non- vanishing components n\ = g^^^^{r), 
"^2 — vdi^- The projections of the Ricci and Riemann tensors on the subspace orthogonal to S 
are 

2r_ 2r+ - 4r_ 

Since i? — for the Reissner-Nordstrom metric, the entanglement entropy of a massless, minimally 
coupled, scalar field {X = 0) and of a conformally coupled scalar field X = ~\R coincide |195] , 

A(I]) 1 .2r+-37_ r+ r 

^RN ^ j^-^ + -rrA )ln— + s(— , 87 

487re^ 90 r+ e r+ 

where A(E) = 47rr^ and s(^) represents the UV finite term. Since s is dimensionless it may 
depend only on the ratio ^ of the parameters which characterize the geometry of the black hole. 

If the black hole geometry is characterized by just one dimensionful parameter, the UV finite 
term in (|87)) becomes an irrelevant constant. Let us consider two cases when this happens. 
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The Schwarzschild black hole. In this case r_ = = 0) and r+ — 2m so that the entropy, 
found by Solodukhin |197] . is 



Historically this was the first time when the subleading logarithmic term in entanglement entropy 
was computed. The leading term in this entropy is the same as in the Rindler space, when the 
actual black hole spacetime is approximated by flat Rindler spacetime. This approximation is 
sometimes argued to be valid in the limit of infinite mass M . We see, however, that, even in this 
limit, there always exists the logarithmic subleading term in the entropy of the black hole that was 
absent in the case of the Rindler horizon. The reason for this difference is purely topological. The 
Euler number of the black hole spacetime is non-zero while it vanishes for the Rindler spacetime, 
respectively the Euler number of the black hole horizon (a sphere) is 2 while it is zero for the 
Rindler horizon (a plane). 

The extreme charged black hole. The extreme geometry is obtained in the limit r_ — >■ r+ 
{q = to). The entropy of the extreme black hole is found to take the form |195j 

Notice that we have omitted the irrelevant constants s(0) and s(l) in (|88|) and (|89|) respectively. 
3.9.2 The dilatonic charged black hole 

The metric of a dilatonic black hole which has mass to, electric charge q and magnetic charge P 
takes the form [120] : 

ds^ = g{r)dT^ + g-\r)dr^ + R^{r)d{de^ + sin^ Odct)'^) (90) 
with the metric functions 

gW- ^'~'^]|;:~''\ R\r)^r^-D\ (91) 
where D is the dilaton charge, D = . The outer and the inner horizons are defined by 



= TO ± Vto2 + D2 - P2 _ ^2 _ (92) 



The entanglement entropy is defined for the outer horizon at r = r^. The Ricci scalar of metric 

dsni) 

p _ on2 (7--r+)(r-r_) 
^- (r2-D2)3 • 

vanishesat the outer horizon, r = r+. Therefore, the entanglement entropy associated with the 
outer horizon is the same for a minimal scalar field {X — 0) and for a conformally coupled scalar 
field (X = -iri?). 



487re^ 90 [r^ — D^) e r+ r+ 



where = 4:Tr{r'^ — D^) is the area of the outer horizon. 
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It is instructive to consider the black hole with only electric charge (the magnetic charge P = 
in this case). This geometry is characterized by two parameters: m and q. In this case one finds 

2 2 2 

r+ = 2m . r_ ~ , rj^ — = 4m(TO ) 

2rn ' 2m 2rn 

so that expression (|93|) takes the form 

SdUaton = + T^(l + 3(1 - -^)) In ^ + s(^) . (94) 

487re^ 180 2m^ e m 

In the extremal limit, 2m^ = g^, the area of the outer horizon vanishes, Ay: — 0, and the whole 
black hole entropy is determined only by the logarithmically divergent ternj^] (using a difi^erent 
"brick wall" method a similar conclusion was reached in |114] ) 

1 r+ 

Sext-dzi = 7^ log — . (95) 
180 e 

In this respect the extreme dilatonic black hole is similar to a two-dimensional black hole. Notice 
that ([95|l is positive as it should be since the entanglement entropy is, by definition, a positive 
quantity. 

The calculation of the entanglement entropy of a static black hole is discussed in the following 
papers [M], [93], [nO], [61], [106], [82], [28], [226], [48], [47], [46], [136], [l35], [175], [M], [l95] . 

[m], [ng, [na, [ni, [ng. 



3.10 Entanglement Entropy of the Kerr-Newman black hole 

The geometry of the rotating black hole is more subtle than that of a static black hole: near the 
horizon the rotating spacetime is no longer a product of a horizon sphere S2 and a two-dimensional 
disk. The other difficulty with applying the technique of the heat kernel to this case is that the 
Euclidean version of the geometry requires the rotation parameter to be complex. Nevertheless 
with some care these difficulties can be overcome and the entanglement entropy of a rotating black 
hole can be computed along the same lines as for a static black hole |170j . In this section we briefiy 
review the results of Mann and Solodukhin |170J . 



3.10.1 Euclidean geometry of Kerr-Newman black hole 

First we describe the Euclidean geometry in the near-horizon limit of the Kerr-Ncwmann black 
hole. The Euclidean Kerr-Newman metric can be written in the form 

dsl = Cdr' + -^^u? + p^de" + sin^ B^^) , (96) 
A — a ) 

where the Euclidean time is t = it and the rotation and charge parameters have also been trans- 
formed a — id, q = iq, so that the metric ()96p is purely real. Here A(r) — (r — f+)(r — f_), where 
f-t = m ± m? + +9^, the quantities uj and Co take the form 

^^SL^±l(dT-dsinHd4>), Co^ ^'' 7" ^ (d0 + , ° „^. dT) (97) 

with fp' = — d^ cos^ 9. This space-time has a pair of orthogonal Killing vectors 

K^dr- " ., ^0 , k ^ dsm^ Odr + (98) 



^Equations l|93l l. II94II and II95II correct some errors in equations (27)-(29) of | 195 | 



26 



which are the respective analogs of the vectors dr and 9^ in the (EucUdean) Schwarzschild case. 
The horizon surface E defined by r = f+ is the stationary surface of the KiUing vector K. Near 
this surface the metric (I96|) is approximately 

dsl = dsl + pldsl^ , (99) 

where p\ — ~ a? cos^ 9 and 

ds| = p\de'^ + -^^-r^ — - sin^ Odi:^ (100) 
P+ 

is the metric on the horizon surface S up to 0{x^), where variable x is defined by the relation 

2 / 

(r — f+) = and 7 = 2^/771^ + + g^. The angle co-ordinate ip — (j) + ° -2-| T and is 

well-defined on I]. The metric ds^^ is that of a two-dimensional disk C2 

dsh, = dx^ + ^dx' (101) 

attached to S at a point {9, ?/;), where x — ''' ^ asin^ </> is an angle co-ordinate on C2. 

Regularity of the metric near the horizon implies the identifications -0 ^ + 27r and x ^ 
X+4:TT'y~^p'^. For this latter condition to hold, independently of 9 on the horizon, it is also necessary 
to identify (r, ip) with (r-|-27r/3i/, (5i—27rJ7 /?//), where ft — ° .3. is the (complex) angular velocity 

and Ph — (f+ — d?)/ ymF+dF+lp . The identified points have the same coordinate ■)/'■ 

Near S we therefore have the following description of the Euclidean Kerr-Newman geometry: 
attached to every point (^,"0) of the horizon is a two-dimensional disk C2 with coordinates (x,x). 
The periodic identification of points on C2 holds independently for different points on the horizon 
S, even though x is not a global coordinate. As in the static case, there is an Abelian isometry 
generated by the Killing vector whose fixed set is E. Locally we have K = d^^. The periodicity 
is in the direction of the vector K and the resulting Euclidean space E is a regular manifold. 

Now consider closing the trajectory of K with an arbitrary period /3 ^ This implies the 
identification (r + 27r/3, (f) — 27rSl/3), and the metric on C2 becomes 

ds^c^ ^=dx^ + a^x^dx^ , (102) 

where x — Pp\{'i% ^ d^)~^x is a new angular coordinate, with period 27r. This is the metric of a 
two dimensional cone with angular deficit 5 — 27r(l — a), a = With this new identification 
the metric (j96p now describes the Euclidean conical space with singular surface E. 

The difference of the Kerr-Newman metric from the static case considered above is that the 
Euclidean space near the bifurcation surface is not a direct product of the surface E and two- 
dimensional cone C2,a- Instead, it is a nontrivial foliation of Ci^a over E. This foliation however 
shares certain common features with the static case. Namely, the invariants constructed from 
quadratic combinations of extrinsic curvature of E vanish E identically. 

3.10.2 Extrinsic curvature of the horizon 

In the case of a static black hole we have argued that the presence of an Abelian isometry with 
horizon being the stationary point of the isometry guarantees that the extrinsic curvature identi- 
cally vanishes on the horizon. In fact this is also true in the case of rotating black hole. The role of 
the Abelian isometry generated by the Killing vector K is less evident in this case. That is why in 
this subsection, following the analysis of ref. |170J . we explicitly evaluate the extrinsic curvature for 
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the Kerr-Newman black hole and demonstrate that quadratic invariants, that can be constructed 
with the help of the extrinsic curvature, vanish on the horizon. 

With respect to the Euclidean metric (j96p we may define a pair of orthonormal vectors {ua = 
, a = 1,2}: 

n[ = J^; n^ = i!i=£l, 4 = ^. (103) 

Covariantly these are 





nl=^-\l—asm''e. (104) 
\ P 

The vectors and are normal to the horizon surface E (defined as r = r_|_, A(r = r_|_) = 0), 
which is a two-dimensional surface with induced metric 7^,^ = — n^nj, — nj^n^- With respect 
to the normal vectors n", a = 1,2 one defines the extrinsic curvatures of the surface E: = 
— 7^7^VQn^. The exact expression for the components of extrinsic curvature is given in [170] . 
The trace of the extrinsic curvature, n"- — K'^^g'^'^, 

= -5j4 , = (105) 
p y p 

vanishes when restricted to the horizon surface S defined by condition A(r = f+) = 0. Moreover, 
the quadratic combinations 

1 f,,_2r^A 2 _ 2a^ cos^ 9 A 

vanish on the horizon E. Consequently, we have kJJ^k"'^" = on the horizon. 
3.10.3 Entropy 

Applying the conical singularity method to calculate the entanglement entropy of a rotating black 
hole we have to verify that i) the curvature singularity at the horizon of a stationary black hole 
behaves in the same way as in the static case and ii) there are no extra surface terms in the heat 
kernel expansion for the rotating black hole. The first point was explicitly checked in |170) : the 
curvature formulas (|56)) - (|59)) are still valid in the stationary case. Regarding the second point, it 
was shown by Dowker [7S] that for a generic metric with conical singularity at some surface E the 
only modification of the surface terms in the heat kernel expansion (|70p are due to the extrinsic 
curvature of E. For example, the surface coefficient may be modified by integrals over E of 
terms k^Kq and kJJj^k"^''. Since, as was shown in the previous section, these terms identically 
vanish for the Kerr-Newman metric there is no modification of the surface terms in this case. The 
expression for the entropy ([82|) thus remains unchanged in the case of rotating black hole. The 
Ricci scalar for the Kerr-Newmann metric is zero, R — 0. The integrals of the projections of Ricci 
and Riemann tensors over horizon surface are 



Rijij = Stt ^ 47r-:2 ln( 



,r 4 



9^ A I (^+ ^ , +a. 
1^ 



Ru = in^ { 1 + ' ln(::^^) ) . (107) 
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The analytic continuation of these expressions back to real values of the parameters a and q requires 
the substitution 



-9 2-^9 2 -X 

q = -q , a = -a , r+ = r+ 

-H^ r =-tan-i — . 108 

a r< — a a r+ 



With these identities the quantum entropy of the Kerr-Newman black hole reads [170] 



5™ . |S) + 1 (i _ + (:l±f!) ta„-.(^„) ,.1 , (109, 

487re^ 45 \ 4r^ ar+ r+ J e 



where A{Yj) — 47r(r^ + a?) is area of the horizon E. In the limit a — > this expression reduces 
to that of the Reissner- Nordstrom black hole (|87)) . An interesting and a still somewhat puzzling 
feature of this result is that, in the case of the Kerr black hole, described by the Kerr-Newman 
metric with vanishing electric charge ((7 = 0), the logarithmic term in the entropy does not depend 
on the rotation parameter a and is the same as in the case of the Schwarzschild black hole. In 
particular for the extreme Kerr black hole (5 = 0, m = a) one has 



5Ke.. = ^ + ^ln^. (110) 

487re2 45 e 



The entropy of the Kerr black hole in the "brick wall" model was calculated in [45] and a result 
different from (|109p was found. The subsequent study in i92. has, however, confirmed ()109p . 



3.11 Entanglement entropy as one-loop quantum correction 

A natural point of view on the entanglement entropy of black hole is that this entropy, as was 
suggested by Callan and Wilczek [33], is the first quantum correction to the Bekenstein-Hawking 
entropjH . Indeed, the Bekenstein-Hawking entropy Sbh can be considered as classical, or tree- 
level, entropy. If we restore the presence of the Planck constant Ti the Bekenstein-Hawking entropy 
Sbh is proportional to 1/h while the entanglement entropy Sent is an Ti^ quantity. The total 
entropy of black hole is then a sum 

5* = Sbh + Sent , (HI) 
where all particles that exist in Nature contribute to the entanglement entropy Sent- 



3.12 The statement on the renormalization of the entropy 

As defined in the previous sections the entanglement entropy is a UV divergent quantity. The 
other well-known quantity which possesses UV divergences is the effective action. The standard 
way to handle the UV divergences in the action is to absorb them into redefinition of the couplings 
which appear in the gravitational action. In four dimensions the gravitational action should also 
include the terms quadratic in the Riemann curvature. The renormalization procedure then is well 
studied and is described in the textbooks (see for instance [l^)- The idea now is that exactly 
the same procedure renormalizes the UV divergences in the entropy. In order to demonstrate 
this statement consider a minimally coupled scalar field. For simplicity suppose that the mass of 

^This statement should be taken with some care. Entanglement entropy is a small correction compared to the 
Bekenstein-Hawking entropy if the UV cutoff 1/e is, for example, of order of few GeV (energy scale of the Standard 
Model). However, the two entropies are of the same order if the cutoff is at the Planck scale. I thank G. 't Hooft 
for his comments on this point. 
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the field vanishes. The bare (tree-level) gravitational action in four dimensions is the sum of the 
Einstein-Hilbert term and all possible combinations quadratic in the Riemann curvature, 



Wgr = j Vgd^x (^-^^1^{R + 2Ab) + c,,bR^ + C2.bRI. + C3.bRI^^p^ , (112) 

where Gb, , As, ci.b, C2,s, cs.b are the bare coupling constants in the gravitational action. 

The UV divergences of the gravitational action are computed by the method of the heat kernel 
using the small s expansion For a minimal massless field {X — in the scalar field equation) 
one finds 

1 f ^ 1 f ^ l/"/l„9 1^2 1d2 



These divergences are removed by standard renormalization of the gravitational couplings in the 
bare gravitational action 

Wgr (Gb , C,,s , As) + Wd^v (e) = Wgr (G Ken) , (114) 

where Gren and Ci^ren are the renormalized couplings expressed in terms of the bare ones and the 
UV parameter e 

111 11 

g;;;-g^ + 12^' '^''■^"-'^■^ + 32^36^"' 

C2,ren-C2.B-^^\ne, C3,.e„ = C3.B + ^ ^ In 6 . (115) 

The tree-level entropy can be obtained by means of the same replica trick, considered in the 
previous sections, upon introduction of the conical singularity with a small angle deficit 27r(l — a), 
S{GB,Ci^B) = (ctda — l)Wgr{a). The conical singularity at the horizon E manifests itself in that 
a part of the Riemann tensor for such a manifold -E^ behaves as a distribution having support on 
the surface E. Using formulas one finds for the tree-level entropy 

S{Gb,c,,b) = -^A(E) - / {Sttci^bR + 4,ttc2,bRu + ^i^c^i^BRijij) ■ (116) 

The Bekenstein-Hawking entropy 5 = ^yl(E) is thus modified due to the presence of i?^-terms 
in the action (|112p . It should be noted that (|116p exactly coincides with the entropy computed 
by the Noether charge method of Wald [220] . [144] (the relation between Wald's method and the 
method of conical singularity is discussed in |140] ). 

The UV divergent part of the entanglement entropy of black hole has been already calculated, 
see (|82p. For a minimal massless scalar one has 

Sd^v = 4^ - f (r- - 2i?,:,„)) Ine . (117) 



487re2 1447r _ 

The main point now is that the sum of the UV divergent part (|117p of the entanglement entropy 
and the tree-level entropy (|116l) 

S{Gb,Ci,b) + Sdivi^) = S{Gr 

takes again the tree-level form (|116p if expressed in terms of the renormalized coupling constants 
Grem Ci,ren defined in (jllSp . Thus, the UV divergences in entanglement entropy can be handled by 
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the standard renormalization of the gravitational couplings. So that no separate renormalization 
procedure for the entropy is required. 

It should be noted that the proof of the renormalization statement is based on a nice property of 
the heat kernel coefficients a„ (j68p on space with conical singularity. Namely, up to (1 — a)^ terms 
the exact coefficient a„ = aj'^'^f on the conical space Ea is equal to the regular volume coefRcient 
aj^"^^ expressed in terms of the complete curvature, regular part plus a delta-like contribution, using 
relations (f55|) 

a„(i?„) = a^^f (i?) + a| = <;s(i? + R'"'^) + 0((1 - af) . (119) 

The terms quadratic in i?'*"'^ are not well defined. However these terms are proportional to (1 — a)^ 
and do not affect the entropy calculation. Thus, neglecting terms of order (1 — a)^ in the calculation 
of entropy the renormalization of entropy (|118p directly follows from the renormalization of the 
effective action (|114p . 

That the leading 1 /e^ divergence in the entropy can be handled by the standard renormalization 
of Newton's constant G has been suggested by Susskind and Uglum [212] and Jacobson (142j . 
That one has to renormalize also the higher curvature couplings in the gravitational action in 
order to remove all divergences in the entropy of the Schwarzschild black hole was suggested by 
Solodukhin [197) . For a generic static black hole the renormalization statement was proved by 
Fursaev and Solodukhin in |111| . In a different approach based on t'Hooft's "brick wall model" 
the renormalization was verified for the Reissner-Nordstrom black hole by Demers, Lafrance and 
Myers [62j . For the rotating black hole described by the Kerr- Newman metric the renormalization 
of the entropy was demonstrated by Mann and Solodukhin [170] . The non-equilibrium aspect (as 
defining the rate in a semiclassical decay of hot flat space by black hole nucleation) of the black 
hole entropy and the renormalization was discussed by Barbon and Emparan ^llj . 



3.13 Renormalization in theories with a modified propagator 

Let us comment briefly on the behavior of the entropy in theories described by a wave operator 
V = F(— V^) which is a function of the standard Laplace operator V^. In flat space this was 
analyzed in section 12.121 As is shown in (183^ there is a precise relation between the small s 
expansion of the heat kernel of operator F(— V^) and that of the Laplace operator — V'^. The 
latter heat kernel has the standard decomposition 

Tre-^ = ^^^^a„."-/^ (120) 

The heat kernel of operator F(— V^) then has the decomposition |183j 

Tre-^(-v^) = ^^5:a„r„(.), (121) 



n=0 



where 



r Pd-2n{s) n<d/2 

Tn{s) = < ^^g^^n-d/2^~sF(g) „ > ^/g ^^^^^ 

I 9=0 

In even dimension d the term 7^/2 (s) — 1. This decomposition is valid both for regular manifolds 
and manifolds with a conical singularity. If a conical singularity is present, the coefficients a„ have 
the standard decomposition into regular a^'^^ and surface parts as in ([551) ■ The surface term for 
n = 1 is just the area of the surface S while the surface terms with n > 2 contain surface integrals 
of {n — l)-th power of the Riemann curvature. Thus ([12ip is a decomposition in powers of the 
curvature of the spacetime. 
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The functions P„ are defined in (|40|) . In particular, if F{q) = {k > 0) one finds that 

= (123) 

The terms with n < d/2 in decomposition (jl2ip produce the UV divergent terms in the effective 
action and entropy. The term n = d/2 gives rise to the logarithmic UV divergence. In d dimensions 
the area term in the entropy is the same as in fiat spacetime (see ea. (l4ip ). In four dimensions 
{d = 4) the UV divergent terms in the entropy are 



5 = 



^ r ^ P2(s) _ J— [ (r^ i(i?^,, _ 2i?„-,„)^ Ine . (124) 
487r J,2 s ^ ' 1447r \ 5^ ^ ^7 ^ ^ 



We note that an additional contribution to the logarithmic term may come from the first term in 
l|124p (for instance, this is so for the Laplace operator modified by the mass term, F{q) = q + m?). 
In the theory with operator i^(— V^) Newton's constant is renormalized as [183] 

' P2[s) (125) 



Gren Gb 127r J ^2 S 

while the higher curvature couplings i — 1,2,3 in the effective action are renormalized in the 
same way as in (jll5l) . The renormalization of G and {ci} then makes finite both the effective 
action and the entropy exactly in the same way as in the case of the Laplace operator — V^. Thus 
the renormalization statement generalizes to the theories with modified wave operator i^(— V^). 

3.14 Area law: generalization to higher spin fields 

In this section we will focus only on the leading UV divergent term, proportional to the area of 
the horizon. The proportionality of the entanglement entropy to the area is known as the "area 
law". As we have discussed already for the case of a scalar field, this term in the entanglement 
entropy of a black hole is the same as in fiat spacetime. In fiat Minkowski spacetime, for a field 
of spin s, massive or massless, including the gauge fields, the calculation of entanglement entropy 
effectively reduces to the scalar field calculation, provided the number of scalar fields is equal to the 
number of physical degrees of freedom of the spin-s field in question. The contribution of fermions 
comes with the weight 1/2. Thus we can immediately write down the general expression for the 
entanglement entropy of a quantum field of spin s in d dimensions, 

where T>s{d) is (with weight 1/2 for fermionic fields) the number of physical (on-shell) degrees of 
freedom of a particle of spin s in d dimensions. For gauge fields this assumes gauge fixing. In 
particular one has 

2[d/2\ 

2?i/2(rf) = (127) 

for Dirac fermions, 

X'i(d) = (d- 2) • dimG, (128) 

for the gauge vector fields (including the contribution of ghosts), where dimG is the dimension of 
the gauge group, 

2[rf/2] 

V^/.M^{d-2)—- (129) 
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for the Rarita-Schwinger particles of spin 3/2 with gauge symmetry (gravitinos) , and 



= (130) 

for massless spin-2 particles (gravitons). 

The result for Dirac fermions was first obtained by Larsen and Wilczek [162] . |163| and later 
in a paper of Kabat jl50] . The contribution of the gauge fields to the entropy was derived by 
Kabat [150j . The entropy of the Rarita-Schwinger spin- 3/ 2 particle and of a massless graviton was 
analyzed by Fursaev and Miel^ |110) . 

3.15 Renormalization of entropy due to fields of diff"erent spin 

The effective action of a field of spin s can be written as 

^tlir ^Tre-^<^' . (131) 
2 J^2 s 

The second order covariant operators acting on the spin-s field can be represented in the following 
general form 

A(") =-V2 + x("', (132) 

where the matrices X^^^ depend on the chosen representation of the quantum field and are linear 
in the Riemann tensor. Here are some examples [13], [U] 

^^R, X%^^ = -^R5ab , = ±R^. , 

(2) 

^fj.u,af:) = -^RidfJ-ag^fi + gfijigva) - Raf^gflu - Rpugafj. - Rnaul3 " Ri^afifi , (133) 

where are gamma- matrices. The coefficient a^*^^ in the small s expansion (|67p - (|69p of the heat 
kernel of operator (|132p has the general form 

a[^^ ^ f {£MR-TrX(^)), (134) 
Je 

where Ds{d) is the dimension of the representation of spin s, 

Z...3/. = ^2[^/^i.....(^^^lM±^. (135) 

Ds (d) can be interpreted as the number of off-shell degrees of freedom. 
Let us consider some particular cases. 



^Among other things the authors of IllOl observe certain, surprising, non-smooth behavior of the heat kernel 
coefficients for the spin-3/2 and spin-2 fields in the limit of vanishing angle deficit. 
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Dirac fermions (s = 1/2). The partition function for Dirac fermions is Z1/2 
In this case TrX(i/2) = i 2[''/2] ^ a,nd hence 

6 jb 

where 2?i/2 was introduced in (|127p . We note that the negative sign in (|136p in combination 
with negative sign for fermions in the effective action (|13ip gives the total positive contribution to 
Newton's constant. The renormahzation of Newton's constant due to Dirac fermions is 

' - ' ' ' ^^/^ ' (137) 



4Gren 4G (4^)^(d-2) 6 

Comparison of this equation with the UV divergence of entropy (|126p for spin- 1/2 shows that 
the leading UV divergence in the entropy of spin-1/2 field is handled by the renormahzation of 
Newton's constant in the same manner as it was for a scalar field. 

The Rarita-Schwinger field (s = 3/2). The partition function, including gauge fixing and the 
Faddeev-Popov ghost contribution, in this case is 

= deti/2A(3/2)det^iA(i/2) ^ (133) 

so that the appropriate heat kernel coefficient is 

ai=af/^'-2a(i/2) = _^ f (139) 

Je 

where I'3/2 is introduced in (|129p . The renormahzation of Newton's constant 

1 _ 1 1 2?3/2 1 



4G,e„ 4G 



(47r)^(d-2) 



(140) 



then, similarly to the case of Dirac fermions, automatically renormalizes the entanglement entropy 
». 

This property, however, does not hold for all fields. The main role in the mismatch between the 
UV divergences in the entanglement entropy and in Newton's constant is played by the non-minimal 
coupling terms X^*) which appear in the field operators (jl32[) . 

3.16 The puzzle of non-minimal coupling 

The simplest case to consider is 

A non-minimally coupled scalar field. In this case one has TrX(°) = ^R, where ^ is the 
parameter of non-minimal coupling. The renormahzation of Newton's constant 



1 1 ,1 1 

(47r)^(d-2)' 



4G.,„ 4G ' (A^\^(d^'?.\^6 ^Kd-^ ^^^^^ 



is modified due to the presence of the non-minimal coupling ^ in the scalar field operator. At 
the same time, the entropy calculation on a Ricci flat background {R = 0) is not affected by the 
non-minimal coupling since the field operator for this background is identical to the minimal one. 
This simple reasoning shows that the area law in the case of a non-minimally coupled scalar field 
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is the same as in the case of the minimal scalar field (|8T|) . Clearly, there is a mismatch between 
the renormalization of Newton's constant and the renormalization of the entanglement entropy. 
One concludes that, in the presence of non- minimal coupling, when the Riemann tensor appears 
explicitly in the action of quantum field, the UV divergence of the entanglement entropy can not be 
handled by the standard renormalization of the Newton'c constant. The mismatch in the entropy 
is 

It is an important fact that there is no known way to give a statistical meaning to this entropy. 
Moreover, (jl42l) does not have a definite sign and may become negative if ^ is positive. In some 
respect this term is similar to the classical Bekenstein-Hawking entropy: both entropies, at least 
in the framework of the conventional field theory, do not have a well-defined statistical meaning. 
There is a hope that in string theory the terms similar to (|142p may acquire a better meaning. 
This question is however still open. 

We should note that on a space with a conical singularity one can consider the Ricci scalar in 
the non-minimal scalar operator as the complete curvature including the (5-like singular term as 
in ([55)) . Then the differential operator — (V^ +£,R) contains a delta-like potential concentrated on 
the horizon surface S. The presence of this potential modifies the surface terms in the heat kernel 
in such a way that [196] 

afiO = «f (e = 0) - 47re(l -a) j^l + 0{l- , (143) 

where of (^ = 0) is the surface term of (|70p without the non-minimal coupling, the term 0(1 — a)^ 
is ill-defined (something like (5^(0)), however it does not affect the entropy calculation. If we now 
apply the replica trick and calculate the entropy corresponding to the theory with the heat kernel 
with the surface term (fT43|) we get that [196], [163], [E] 

Sd^v = - (144) 

(d- 2)(47r) — £^^-2 6 

This divergence takes the form consistent with the UV divergence of Newton'constant (|141l) . How- 
ever, we can not interpret this entropy as a contribution to the entanglement entropy since the 
presence of the delta-like potential in the Euclidean field operator is not motivated from the point 
of view of the original Lorentzian theory, for which the entanglement entropy is calculated. More- 
over, (I144p is not positive if ^ > 1/6 while the entanglement entropy is supposed to be a positive 
quantity. 

Similar features are shared by other non-minimally coupled fields. 



Abelian vector field. After gauge fixing the partition function of an Abelian gauge field is 

Z = det-i/2^^i) .^^^.^(0)^ (^45) 

where A^^-* is operator defined in (|132p with sign -I- in the matrix xj^J (|133p . For the effective 
action Weff = — In Z we find. 




(146) 
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where 

(s=i) (s=o) T\ f T} nA'7\ 

ai = al ' ~2al ' = ( — 1) J R, (147) 

and 'Di{d) ~ d — 2 is the number of on-shell degrees of freedom of the Abelian vector field. The 
renormalization of Newton's constant is 

' - ' ' ' ("8) 



4G 

ren 

4G (4^)V(d_2)^ 6 



Comparison with (|126l) shows that there is again a mismatch between the UV divergences in the 
entropy and in Newton's constant. This mismatch originates from the non- minimal term xj^J in 
the Laplace type field operator for the vector field. 



Massless graviton. The partition function of a massless graviton in d dimensions, after gauge 
fixing and adding the Faddeev-Popov ghost contribution, is 

Z = det-i/2A(2) . det aL'^ • det-^/^ACo) ^ (^49) 

where A^^ is operator defined in (|132p with sign — in the matrix xj^J (|133p . The operator A'^^^ 
governs the dynamics of the tensor perturbations which satisfy the condition ^^{h^u — 29 

h) = 0, 

the operator A^^-* is due to the contribution of the conformal mode while the determinant of 
operator aL^-* is due to the Faddeeev-Popov ghosts. Hence one has in this case that 

^.=2 o ^ _ Md) f ^ , _ d^-d-' ^ 



ai = - 2a\^-^> + a\^-'" = (^-^ - c(d)) R , c{d) = ^— , (150) 

where I?2(d) = ^id-s) ^j^^ number of on-shell degrees of freedom of a massless spin-2 particle. 
The renormalization of Newton's constant is 

4:Gren 4G {An)— {d - 2) 6 ^ 

Again, we observe the mismatch between the UV divergent terms in the entropy (|126p and in 
Newton's constant, this time due to the graviton. 

To summarize, the UV divergences in the entanglement entropy of minimally coupled scalars 
and fermions are properly renormalized by the redefinition of Newton's constant. It happens that 
each minimally coupled field (no matter bosonic or fermionic) contributes positively to Newton's 
constant and positively to the entropy of black hole. The contributions to both quantities come 
proportionally that allows the simultaneous renormalization of both quantities. The mismatch 
between the UV divergences in the entropy and in Newton's constant appears for gauge bosons: 
the Abelian (and non- Abelian) vector fields and gravitons. The source of the mismatch are those 
non-minimal terms A'^'*) in the field operator which contribute negatively to Newton's constant 
and do not make any contribution to the entanglement entropy of the black hole. At the time this 
review is being written, the appropriate treatment of the entropy of non-minimally coupled fields 
is not yet available. 

3.17 Comments on the entropy of interacting fields 

So far we have considered free fields in a fixed gravitational (black hole) background. The in- 
teraction can be included by adding a potential term J d'^xV{'ip) to the classical action, here 
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ip — {-0* J i — 1,..,A^} is set of fields in question. In the one-loop approximation one splits 
= V'c + V'g : where V'c is the classical background field and ■i/'g is the quantum field. The integra- 
tion over '0g then reduces to calculation of the functional determinant of operator 2? -I- M^(?/'c), 
where M^?- = dfjV{tlJc)- The fields i/'c representing the classical background are in general func- 
tions on the curved space-time. In some cases these fields are constants that minimize the potential 
V{ip). The matrix Mij{x) plays the role of an a;-dependent mass matrix. In the approximation 
when one can neglect the derivatives of the matrix M the heat kernel of operator V + is 
presented as the product Tre^''^ • Tre^*^^ . Using the already calculated trace of the heat kernel 
rpj-g-sP Qj-^ space with a conical singularity one obtains at one-loop the entanglement entropy of 
the interacting fields. In d dimensions one obtains j204j 

1 f d 



12(47r)— Jt. 2 

where we used that /J^ s~''/^e^*^ — T{1 — |, M^e^). Using the asymptotic behavior r(— a, a;) = 
a~^a;~" -|- .., we find that the leading UV divergence of the entropy (|152[) is again (multiplied by 
N) (|8T|) and is thus not affected by the presence of the interaction in the action. The interaction 
however shows up in the sub-leading UV divergent and the UV finite terms. For instance, in four 
dimensions on a flat background we find 

^ = + 4^X ((7- 1)™^ • (153) 

We see that the leading UV divergent term proportional to the area is not modified by the presence 
of the self- interaction. The mass matrix is function of the background field ijic- The result 
(|153p thus indicates that at tree-level the entropy should contain terms additional to those of the 
standard area law, which depend on the value of field "0 at the horizon. In order to illustrate this 
point consider a model of a single field (a two-dimensional model of this type was considered 
in [151J . in four dimensions the role of self-interaction was discussed in [46] ) 



d^xy^ (^(V0)2 + + 



(154) 



^-^(«-01ne, (155) 



where we included the term with the non-minimal coupling. In fact if we had not included this 
term, it would have been generated by the quantum corrections due to the self-interaction of the 
field V'^- This is a well-known fact, established in f l81| . The renormalized non-minimal coupling 
in the model (|154p is 

-(- 

where we omit the terms of higher order in A. Splitting in (|154p the field tp into classical and 
quantum parts we find AP — + Xip"^. Suppose for simplicity that the background metric is flat. 
Then to leading order in A the entanglement entropy (omitting the UV finite terms) is 

5.„ = ^ + ^£v.?ln.. (156) 

This entropy should be considered as a quantum correction to the tree-level entropy 

Stree = - ^< f , (157) 

which follows from the action Wgr -\- W[il>\. We see that the logarithmic divergences in (|156p 
and (|155l) agree if ^ = 0. On the other hand, the renormalization of ^ (|155p does not make the 
total entropy Stree + Sdiv completely UV finite. This is yet another manifestation of the puzzling 
behavior of the non-minimal coupling. 
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4 Other related methods 



4.1 Euclidean path integral and thermodynamic entropy 

In 1977 Gibbons and Hawking [119] developed a method based on the Euchdean path integral for 
studying the thermodynamics of black holes. In this method one obtains what may be called a 
thermodynamical entropy. One deals with metrics which satisfy the gravitational field equations 
and thus avoids the appearance of metrics with conical singularities. The entanglement entropy, on 
the other hand, has a well-defined statistical meaning. In ordinary systems the thermodynamical 
entropy and the statistical (microscopical) entropies coincide. For black holes the exact relation 
between the two entropies can be seen from the following reasoning [196\ 

Consider a gravitationally coupled system (gravity plus quantum matter fields) at some arbi- 
trary temperature T = (/3)~^. A standard way to describe a thermal state of a field system is to 
use an Euclidean path integral over all fields in question defined on manifold with periodicity 27r/3 
along the time-like Killing vector. Suppose that it is a priori known that the system includes a 
black hole. Thus there exists a surface E (horizon) which is a fixed point of the isometry generated 
by the killing vector. This imposes an extra condition on the possible class of metrics in the path 
integral. The other condition to be imposed on metrics in the path integral is the asymptotic be- 
havior at infinity: provided the mass M and the electric charge Q of the gravitational configuration 
are fixed, one has to specify the fall-off of the metrics for large values of r. Thus, the Euclidean 
path integral is 

Z{P, M,Q) = J Vg^, j 2?Ve-^-[9l+^"-[^'3] , (158) 

where the integral is taken over /3-periodic fields ^P{t,x^) = ip{T + P,x) and over metrics which 
satisfy the following conditions: 

i) g^^^ possesses an Abelian isometry with respect to the Killing vector dr] 

ii) there exists a surface S (horizon) where the Killing vector dr becomes null; 

iii) asymptotic fall-off of metric g^u at large values of radial coordinate r is fixed by the mass M 
and electric charge Q of the configuration. 

Since the inverse temperature /3 and mass M in the path integral are two independent param- 
eters, the path integral ()158|) is mostly over metrics which have a conical singularity at the surface 
E. The integration in ()158|) can be done in two steps. First, one computes the integral over matter 
fields ■0 on the background of a metric which satisfies conditions i), ii) and iii). The result of this 
integration is the quantity ([T5|) used in the computation of the entanglement entropy, 

Semiclassically, the functional integration over metrics in (|158p can be performed in a saddle-point 
approximation, 

Z(/3, Q, M) = e-i^""['''ff('5)l , (160) 
where metric (/3) is a solution to the saddle-point equation 

^-}^^ldM = Q^Wtot = Wgr[l3,g] + W[l3,g], (161) 



^As I have learned recently (private communication of R. Myers), Jacobson and Myers (unpublished) had similar 
ideas back in the 90-s. 
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whith the inverse temperature /3 kept fixed. The solution of this equation is a regular (without 
conical singularities) metric g^v{(3). This is an on-shell metric which incorporates the quantum 
corrections due to the vacuum polarization by the matter fields. It can be also called an equilibrium 
configuration which corresponds to the fixed temperature /3~^. In the saddle point approximation 
there is a constraint relating the charges at infinity M and Q and the inverse temperature /3: 
/3 = /3(M,g). 

The thermodynamic entropy is defined by the total response of the free energy F — —(3 ^ In Z{(3) 
to a small change of the temperature, 

Std = iS^dpF = (/3d/3 - l)WtotW,gp) (162) 

and involves, in particular, the derivative of the equilibrium configuration g^^iP) with respect to 

apWtot = 0|3Wtot[p,g\^ 7 ■ 163) 

For an equilibrium configuration, satisfying (|16ip . the second term in (|163p vanishes and thus the 
total derivative with respect to /3 coincides with a partial derivative. 

Thus, in order to compute the thermodynamical entropy one may proceed in two steps. First, 
for a generic metric which satisfies the conditions i), ii) and iii) compute the off-shell entropy using 
the replica method, i.e. by introducing a small conical singularity at horizon. This computation is 
done by taking a partial derivative with respect to (3. Second, consider this off-shell entropy for an 
equilibrium configuration which solves equation (|16ip . Since for the classical gravitational action 
(I112p one finds {(3dp — l)Wgr[P,g] = S{GB,CiB) (|116p and for the quantum effective action one 
obtains the entanglement entropy (/39/3 — l)V[^[/3, = Sent the relation between the entanglement 
entropy and thermodynamical entropy is given by 

Std = S{Gb,c.,b) + Sent ■ (164) 

Therefore, the entanglement entropy constitutes only a (quantum) part of the thermodynamical 
entropy of the black hole. The thermodynamical entropy is defined for equilibrium configurations 
satisfying the quantum corrected Einstein equations (jl61l) . Thus these configurations are not 
classical solutions to the Einstein equations but incorporate the quantum (one-loop) corrections. 
These configurations are regular metrics without conical singularities. The UV divergences in 
the free energy for these configurations are renormalized in a standard way and thus for the 
thermodynamical entropy the renormalization statement discussed above holds automaticall}|f|. 

In flat spacetime the quantum (one-loop) thermodynamical and statistical entropies coincide 
as was shown by Allen [5] due to the fact that the corresponding partition functions differ by terms 
proportional to (3. In the presence of black holes the exact relation between the two entropies has 
been a subject of some debate (see for example [86], |198j ). The analysis made in |106j however 
shows that in the presence of black hole the Euclidean and statistical free energies coincide provided 
an appropriate method of regularization is used to regularize both quantities. 



4.2 't Hooft's brick wall model 

In 1985 t' Hooft ^213j proposed a model which was one of the first successful demonstrations that 
an entropy that scales as an area can be associated, in a rather natural way, with a black hole 
horizon. The idea of t' Hooft's calculation was to consider a thermal gas of Hawking particles 

^This is true for minimally coupled matter fields. In the presence of non-minimal couplings there appear extra 
terms in the thermodynamical entropy which are absent in the entanglement entropy, as we discussed earlier. 
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propagating just outside the black hole horizon. The entropy in the canonical description of the 
system is calculated by means of the WKB approximation. Provided the temperature of the gas 
is equal to the Hawking one the result of this calculation is unambiguous. There is however an 
important subtlety: the density of states of a Hawking particle becomes infinite as one gets closer 
to horizon. The reason for this is simple. Close to the horizon all particles effectively propagate in 
the so-called optical metric. The later is conformally related to the black hole metric 

dslH = ''9ir)df + g~\r)dr^ + r^duJd-2 (165) 

as follows 

dslpt = 9'\r)d.3lH = -dt^ + 9'\r)dr' + r'g-\r)dul_^ , (166) 

where dw^_2 is the metric of the {d — 2)-unit sphere. In the optical metric the near-horizon region, 
where the metric function in (|165p can be approximated as 

PH 

occupies an infinite volume. Clearly, the infinite volume contains an infinite number of states. In 
order to regularize this infinity t' Hooft introduced a brick wall, an imaginary boundary at some 
small distance e from the actual horizon. The regularized optical volume then is divergent when e 
is taken to zero 

Vopt^n,^2 f drr^-^g-"/^ ^ A{^)^, (167) 

where A{Yi) = r'!i^^Vld^2 is the area of the horizon and e ^ y^/ffTC^^T^^^^y is the invariant distance 
between the brick wall {r = r^) and the actual horizon (r = r^^). The entropy of a gas of massless 
particles at temperature T = [3^^ confined in the volume Vopt in d spacetime dimensions 

Sbw ~ Vopt^^^'* ^ (168) 

in the optical metric, is proportional to the horizon area. We should note that the universal 
behavior of the regularized optical volume (|167p in the limit of small e and its proportionality to 
the horizon area in this limit was important in establishing the result ()168|) . 



4.2.1 WKB approximation, Pauli-Villars fields 

In the original calculation by 't Hooft one considers a minimally coupled scalar field which satisfies 
the Klein-Gordon equation 

(V2-to2)(^ = (169) 

on the background of a black hole metric (jl65p and imposes a brick wall boundary condition 

ip{x) = at r = re . (170) 

Consider for simplicity the four dimensional case. Expanding the scalar field in spherical coordi- 
nates tf = e~'^'^*Yi,rn{0,4')f{r) one finds that equation (|169|) becomes 

uj'g-\r)fir) + r'^ dr{r^ g{r)dr f {r)) - (^^ + m')f{r) = . (171) 
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One uses the WKB approximation in order to find a solution of this equation. In this approximation 
one represents f{r) = p{r)e^^''^\ where p{r) is a slowly varying function of r while S{r) is a rapidly 
varying phase. One neglects derivatives of p(r) and the second derivative of S{r) on obtains the 
radial function in the form 



/(r) =p(r)e±*r#)'=('''''^\ k{r, l,E) = ^E^- {m^ + ^l±ll)g{r) , (172) 

valid in the region where fc^(r) > 0. The latter condition defines a maximal radius r^^^i which is a 
solution to the equation k'^{r^j) ~ 0. For a fixed value of the energy E, by increasing the mass m 
of the particle or the angular momentum /, the radius r^^,; approaches r+ so that the characteristic 
region where the solution (jl72l) is valid is in fact the near horizon region. One imposes an extra 
Dirichlet condition tp = ai r — r^^^i so that the one-particle spectrum becomes discrete 

—-k{r,u},l)=27rn, (173) 

where n is an integer. This relation is used to count the number of one-particle states that 
correspond to fixed values of energy w and angular momentum I, 

1 r^-' dr 

n{LO,l)^- k{r,u:,l). (174) 

Calculating the total number of states which have the same energy E, one has to sum over I. This 
sum can be approximated by an integral 

n{uj)^jdl{2l + l)n{u:,l)^'^jj^k\r,u:), k{r,u:) = ^Juj^ - m^g{r) , (175) 

where r^^ is determined by condition that k{r,uj) = 0. 

In the near horizon region one approximates the metric function g{r) in (|166p by the first two 
terms in the expansion in powers of (r — r_|_), 

Air 

9{r)^—{r~r+) + C{r~r+)\ (176) 

PH 

where is the inverse Hawking temperature and r+ is the horizon radius. Constant C is related to 

2 

the curvature of spacetime near the horizon. The radial position of the brick wall is = r-|_ -I- , 
where e is the geodesic distance between the brick wall and the horizon. Focusing only on the brick 
wall divergent terms, one obtains for the number of states (|175|) 

^ l4^ + 1^4^ (^^^ - ^) + j ^ ■ ^'''^ 

In a thermal ensemble of scalar particles at fixed temperature T = /3^^, each state in the 
one-particle spectrum can be occupied by any integer number of quanta. One gets for the free 
energy 

/3F= r doj'^^Hl-e-f'n (178) 

Jr. dbJ 



or, integrating by parts, 



PF = -pl -^d^. (179) 
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Substituting here ()177|) and using the integrals 



e/3<^ - 1 6/32 ' 7q eP'^ - 1 15/34 
one calculates the divergent terms in the free energy 



(180) 



360e2 /34 ^360"^ r+'/J" 12 
and, using equation S — 0^dpF\^=pjj, the entropy 

Due to the relations 

Crl - ^ - ^ / (^.. - 2%»j) and 47rr^ = / 1 (183) 
this expression for the entropy can be rewritten in a completely geometric form 

The leading term proportional to the area was first calculated in the seminal paper of 't Hooft 
[Hg. The area law in the brick wall model was also studied in [T73], [HI], 0, [IH], [TTT] . 

It is an important observation made by Demers, Lafrance and Myers in [62] that the brick wall 
divergences are in fact the UV divergences. This can be seen in the Pauli-Villars regularization as 
was first done in |62) . Applying the Pauli-Villars regularization scheme for the four-dimensional 
scalar field theory studied here, one introduces five regulator fields {ipi, i = 1,...,5} of different 
statistics and masses {m^, i = 1,...,5} dependent on the UV cut-off /x [52]. Together with the 
original scalar ipo = ip (mo = m) these fields satisfy two constraints 

5 5 
^ Ai = m and ^ A,m^ = , (185) 

where = -1-1 for the commuting fields, and A, = —1 for the anticomuting fields. Not deriving 
the exact expressions for m^, we just quote here the following asymptotic behavior 

5 2 

EA 2i 2 2i 2i 2 J 

AiTTij Inm^ — bi + m In — + m 02 , 

5 2 

VAanm^^ln^ , (186) 

(where bi and 62 are some constants), valid in the limit /i — )■ 00. The total free energy is the sum 
of all contributions, from the original scalar field and the regulators 

5 

/3F=/3^A,FV (187) 

i=0 

It is clear that due to the constraints (|185p all brick wall divergences (with respect to the parameter 
e) in the free energy (|187p and in the entropy cancel. On the other hand, both the free energy and 
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the entropy become divergent if the Pauh-Villars regulator /i is taken to infinity, thus confirming 
their identification as UV divergences. For the free energy one finds 

F = --^riy ^,M^\nM^ - -^^rlcy l^dnM^ . (188) 
24 «2 + zl^ ' » * 1440 /34 + ^ * ' ^ ' 

' i=0 ^ i=0 

and for the entropy one has 

5 5 

i=0 i=0 

Several remarks are in order. 

1) Comparing the entropy calculated in the brick wall model, (|184p or (|189p . with the en- 
tanglement entropy (|82p we see that the structure of the UV divergent terms in two entropies is 
similar. The logarithmic terms in (|82|) and (|189p (or (|184p ) are identical if the black hole metric 
has vanishing Ricci scalar, R ~ 0. This is the case, for example, for the Reissner-Nordstrom black 
hole considered in f62| . The logarithmic terms in the two calculations are however different if the 
Ricci scalar is non-zero. This discrepancy appears to arise due to certain limitations of the WKB 
approximation. In the exact solutions, known explicitly, for example, for a scalar field in a constant 
curvature spacetime, the mass m always appears in combination — ^R. This, however, is not 
seen in the WKB approximation (jl72p . In fact, if one makes this substitution everywhere in the 
above brick wall calculation, the Ricci scalar would appear in the brick wall entropy in a manner 
which agrees with the entanglement calculation (|82|) . Moreover, an alternative calculation [106 of 
the density of states which does not make use of the WKB approximation results in an expression 
for the entropy which agrees with (|82|) . 

2) The similarity between the two entropies suggests that the UV divergences in the brick wall 
entropy (jl89p can be renormalized by the renormalization of the couplings in the gravitational 
action in the same way as for the entanglement entropy. That this indeed works was demonstrated 
in [62]. 

3) For a non-minimally coupled field we have the same problem as in the case of the entangle- 
ment entropy. For metrics with R = the non-minimal coupling does not show up in the scalar 
field equation and does not change the density of states, the free energy and entropy. On the 
other hand, the non-minimal coupling affects the renormalization of Newton's constant, even if the 
background metric is Ricci fiat. An attempt was made in |199) to modify the Dirichlet boundary 
condition at the brick wall and replace it by a more sophisticated condition, which would depend 
on the value of the non- minimal coupling ^, so that the resulting entropy would have the UV 
divergences consistent with the renormalization of Newton's constant. This attempt however can 
not be considered as successful since it does not reproduce the expected behavior of the entropy 
for large positive values of ^. 

The calculation of the brick wall entropy for a rotating black hole is more complicated due to 
the presence of the superradiance modes in the spectrum. This issue is considered in papers [45] . 
[92], [154], m, [US], [223], [132], [MS], [153]. 

4.2.2 Euclidean path integral approach in terms of optical metric 

Field equation in optical metric. Consider a slightly more general equation than ()169p . by 
including a non-minimal coupling, 

i-V'^ +CR + m'^)^ = (190) 

in the background of the black hole metric g^i/, which takes the form (|165p . The optical metric 
is conformally related to the black hole metric, gfj,i, = e^'^ g^iv, where e^"^ = l/|(7tt| (in the metric 
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(|165p we have that gu = 5 (?"))■ The equation ()190|) can be rewritten entirely in terms of the optical 
metric g^v (|166p as follows 

V = ( i^R + ni^) - i^^V^a - ^-^-^(ya)A , (191) 



where tpopt — e i "^ip, Aopt is the Laplace operator for spatial part 7°^"* = gij of the optical 
metric while the scalar curvature R and the covariant derivative V are defined with respect to 
original metric = e~^'^5^i,. We notice that, since e~^'^ = g{r), the effective potential V in (|19ip 
vanishes at the horizon. This is a general feature of wave equations in the black hole background: 
the fields become, effectively, massless in the near horizon region. The frequency lo which appears 
in equation (|17ip in the brick wall calculation is thus an eigenvalue of the operator H^, 

H^^^ = L^V- ■ (192) 
The canonical free energy and Euclidean path integral. The canonical free energy (I178P 

F = /3-i^n(w)ln(l-e-''"), (193) 

UJ 

where lo are eigenvalues of the spatial operator , n{oj) is the degeneracy of the energy level uj, 
can be represented in terms of the Euclidean path integral for a field theory with wave operator 
{d^ + H^), provided that the Euclidean time r is a circle with period (3. (This property was first 
clearly formulated by Allen [2].) It order to see this in a rather elementary way, we first notice 
that 



ln(l - e-^n ^^B^+lnPu: + J2 + ^) • (194) 

fc=l 



The sum in this expression can be rewritten as a difference of two sums 

fc=l fc=-oo ^ fc=l '"^ 

Each of these sums should be understood in terms of the zeta- function regularization. In particular, 
using the properties of the Riemann ^-function, we find 

°o 4-21,2 J 00 Zl7r2l.2 

fc=l ^ k=l ^ 

Collecting together (fTM)) . ([T95|) and (fT96|) one obtains that 

^^ = -^E"H- + ^E"H E M^^ + ^) (197) 

U) LO k = ~oo 

The second term in (|197p can be expressed in terms of the Euclidean path integral. This can be 
seen as follows. In the Euclidean formulation one first makes a Wick rotation of time t — > —ir. 
The effective action Wopt then is defined by means of the Euclidean path integral 



e 
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At finite temperature one closes the Euclidean time by identifying r and r + /3. The eigen values of 
operator then are i^k, where k - 
of the logarithm of the determinant 



operator then are i^k, where k = 0, ±1, ±2, ... The effective action can be expressed in terms 



1 1 Air'^k^ 

w /c= — oo 

Comparing this with (|197p and defining the vacuum energy as 

UJ 

one arrives at an expression for the free energy f|197p 

F = r^Wopt - Eo . (201) 

Evaluation of the effective action in the optical metric. The effective action Wopt (|198l) 
can be calculated using the heat kernel method. One notes that the heat kernel of the operator 
—d^ + H^ takes the form of a product of two heat kernels, for commuting operators and H^. The 
heat kernel for operator is computed explicitly, provided the periodicity condition, r — > r + /3, 
is imposed. One finds for the trace 

q2 B fc2^2 

Tre"^- = Y . (202) 

ft! — — OO 

So that the effective action takes the form 

^ fc— — oo 

The operator — —Aopt + V (|19ip is defined for the {d — l)-dimensional metric 7j°f*, the spatial 

part of the optical metric (|166p . The trace of the heat kernel of operator = —Aopt + V (|19ip 
can be represented as a series expansion in powers of s, 

TrKfj, = (^J + s J V^tilRopt - V) + 0{s^)^ , (204) 



{■ins) 

where the integration is taken over the spatial part of the optical metric ()166|) . Ropt is the Ricci 
scalar of {d — l)-metric 7°''*. 

For n ^ the integration over the proper time s in (j203l) is regularized for small s due to the 
thermal exponential factor, so that the UV regulator e can be removed. Interchanging the sum 
and the integral one obtains 

Y ^.™-fe-^ = (f)2™A(d-2m)r(^-m), m = 0, 1, 2, .. (205) 



n=l 



Only the term with n = in (|203p contains the UV divergences. This term in the effective action 
is proportional to the inverse temperature /3 and thus it does not make any contribution to the 
entropy. On the other hand, the n = term gives the free energy at zero temperature. Thus, only 
the zero temperature contribution to the free energy is UV divergent as was shown by Dowker and 
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Kennedy |72j . In fact the usual way to renormalize the free energy is to subtract the n — term 
in (|203p or, equivalently, to subtract the zero temperature free energy, — F{T) — F{T — 0). 
With this regularization and using (I205P one obtains a sort of high temperature expansion of the 
effective action. In d dimensions one finds for the regularized action [73] . |74j 

^opt = -^ad)T(^)Vo,t - ^C(d - 2)r(^ -l)j ^t{\Ro,t ~ V) + 0{P'-') , (206) 

where for d = 4 the term 0(/3'*~'') contains also logarithmic term ln/3. In four dimensions the 
first two terms in (|206p are the only terms in the effective action which are divergent when the 
integration in the optical metric is taken up to the horizon. The first term in (I206P and the 
respective term in the free energy and entropy is the contribution of a thermal gas in {d — 1) 
spatial volume Vopt at temperature T = in flat spacetime. The other terms in (|206|) are 
curvature corrections to the flat spacetime result as discussed by Dowker and Schofleld [73], [7i] . 

Let us focus on the four-dimensional case. Defining the regularized free energy F^ = /3~^W^( 
and entropy Sopt = fi'^dpF^ one finds (provided one imposes the condition (3 — fi^, after taking 
the derivative with respect to /3) 

Sl% = "^TliVopt + ^Th J V7^(ii?opt - V) , (207) 

where we omit terms which are finite when the volume integration is extended to the horizon. The 
important observation now is that 

/ V^ilRopt -^) = J V^te-^"iil - OR - ™') : (208) 

where e^^"" = g{r) and R is the scalar curvature of the original black hole metric. We recall that 
the latter is conformally related to the optical metric, 17°^* = e^°'(7^,y, so the relation between the 
scalar curvature in two spacetimes is 

R^pt ^ e-^'^{R - 6VV - 6(Vcr)^) . (209) 

Using this relation and the form of the potential term V one arrives at (|208p . 

2 

Introducing the cut-off e as before, — r+ + one finds for the volume in the optical metric 



+ 4ri I (-R" - 2i?,,,, ) Ine + 0(e) (210) 



and 



j ^t{\Ropt-V = {{\-0R{r+)-m^)I^A+\ne-^ + O{^, (211) 

where = 47rr^ is the horizon area and R{r+) is the value of the scalar curvature at the horizon. 
Putting everything together, one finds that the entropy in the optical metric is 

= 365^? + X {ik^^- - ^^-^^^ -^A" - -^0 1- ■ 

Comparing this result with the entanglement entropy computed earlier, we find complete 

agreement. Notice that e here is in fact the IR regulator, the brick wall cut-off, which regularizes 
the integration in the radial direction in the optical metric. As in t' Hooft's original calculation. 
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this divergence can be transformed into a UV divergence by introducing the Pauli-ViUars regulator 
fields of characteristic mass The UV divergences of the entropy when ^ is taken to infinity then 
are of the same type as in (|212p . This was analyzed in [106' . 

The Euclidean path integral approach in the optical metric was considered in |16j . [58] . [57] . 

i, [He]: W\- 

We remind the reader that the entanglement entropy ([82|) is obtained by using the Euclidean 
path integral in the original black hole metric with a conical singularity at the horizon (for (3 ^ 13 h)- 
The black hole metric and the optical metric are related by a conformal transformation. This trans- 
formation is singular at the horizon and in fact produces a topology change: there appears a new 
boundary at r = r-(- in the optical metric which was a tip of the cone in the original black hole 
metric. Because of this singular behavior of the conformal transformation, the exact relation be- 
tween the two Euclidean path integrals is more subtle than for a regular conformal transformation. 
That the UV divergences in the entropy calculated in these two approaches coincide suggests that 
the equivalence between the two approaches might extend to the UV finite terms. Although in 
arbitrary dimension this equivalence may be difficult to prove, the analysis in two dimensions |198j 
shows that the entanglement entropy and the brick wall entropy are indeed equivalent. This is, of 
course, consistent with the formal proof outlined in Section [3.61 

5 Some particular cases 

In this section we shall consider some particular examples in which the entanglement entropy, 
including the UV finite terms, can be calculated explicitly. 

5.1 Entropy of 2d black hole 

In two dimensions the conformal symmetry plays a special role. This has many manifestations. In 
particular, the conformal symmetry can be used in order to completely reproduce, for a conformal 
field theory (CFT), the UV finite part of the corresponding gravitational effective action. This is 
done by integration of the conformal anomaly. For regular two-dimensional spacetimes the result 
is the well-known non-local Polyakov action. In the presence of a conical singularity the derivation 
is essentially the same although one has to take into account the contribution of the singularity. 
Consider a two-dimensional CFT characterized by a central charge c. For a regular two-dimensional 
manifold the Polyakov action can be written in the form 

WplIM] = ^ //^(^^)' + V^^) ' (213) 

where the field equation for the field ip is V^-i/' — R- On a manifold M" with a conical singularity 
with angle deficit 5 = 2n{l — a) the Polyakov action is modified by the contribution from the 
singularity at the horizon S (which is just a point in two dimensions) so that [197] . |111| 

M^pl[M"] = WplW/Y] + ^(1 - a)^h + 0(1 - af , (214) 

where tph is the value of the field -0 on the horizon. Applying the replica method to the Polyakov 
action ([214p one obtains that the corresponding contribution to the entanglement entropy from 
the UV finite term in the effective action is 

Sftn = J^Tph ■ (215) 

This result agrees with a derivation of Myers |178] who used the Noether charge method of Wald 
|220l in order to calculate the entropy. The easiest way to compute the function ijj is to use the 
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conformal gauge g^i, — e^" 5^^ in which ijj — 2a. Together with the UV divergent part, the complete 
entanglement entropy in two dimensions is 

S^^<^h + l\^-, (216) 
6 6 e 

where A is an IR cut-off. 

Let the black hole geometry be described by a 2d metric 

dsl^ = f{x)dT^ + -^^dx^ (217) 

where the metric function f{x) has a simple zero at a; = x+. Assume that this black hole is placed 
inside a box of finite size L so that x+ < x < L. In order to get a regular space one closes the 
Euclidean time t with period jin , Ph — f'\x+) ■ '^^ easy to see that (|217l) is conformal to the flat 

dx 



disk of radius zq (Inz = ^ TT^-' 



dsl,^^ e^^z^idz^ + z^df"^) , (218) 
1 , ^ 2tt r dx Ph 



2 ' PhJl fix) 2^zo' 

where f = (0 < f < 27r), < z < 1. So that the entanglement entropy of the 2d black hole 
takes the form [HH], [SS] 

12 }{x) 6 e 

where we omit the irrelevant term that is a function of (A, zq) but not of the parameters of the 
black hole and have retained dependence on the UV regulator e. 

As was shown in |198j , the entanglement entropy (j219p is identical to the entropy of the thermal 
atmosphere of quantum excitations outside horizon in the "brick wall" approach of 't Hooft [213] . 

The black hole resides inside a finite size box and L is the coordinate of the boundary of the 
box. The coordinate invariant size of the subsystem complimentary to the black hole is iinv = 

dxj ^ f{x). Two limiting cases are of interest. In the first, the size of the system Li^v is taken 
to infinity. Then, assuming that the black hole space-time is asymptotically fiat, we obtain that 
the entanglement entropy (|219|) approaches the entropy of the thermal gas, 

CTT 

S = y ilnvTff . (220) 

This calculation illustrates an important feature of the entanglement entropy of a black hole placed 
in a box of volume V. Namely, the entanglement entropy contains a contribution of the thermal 
gas that, in the limit of large volume in dimension d, takes the form ([7]). This is consistent with the 
thermal nature of the reduced density matrix obtained from the Hartle-Hawking state by tracing 
over modes inside the horizon. 

The other interesting case is when iinv is small. In this case we find the universal behavior 



5 ^ _ l^ln + :^Li^^ + 0{Li^^)j . (221) 

The universality of this formula lies in the fact that it does not depend on any characteristics of 
the black hole (mass, temperature) other than the value of the curvature R{x+) at the horizon. 
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Consider two particular examples. 

2 

2d de Sitter spacetime is characterized by the metric function f{x) — I — jr, the Hawking 
temperature Th = 1/2ttI. In this spacetime the size of the box is bounded from above, Linv < ttL 
The corresponding entanglement entropy 

SdS2 = ^ In ( tan(rjf^L„„) ) (222) 

6 \nlHe J 

is a periodic function of Linv 

The string inspired black hole [221], |169| is described by the metric function f{x) = 1 — 
e"'*'^. It described an asymptotically flat spacetime. The Hawking temperature is Th = j^. The 
entanglement entropy in this case is 

Sstr = 7 In f 7r4^ smh{2TrTHL,nv)) ■ (223) 

The entropy in these two examples resembles the entanglement entropy in flat spacetime at 
zero temperature (|11|) and at a finite temperature ([T2|) respectively. 



5.2 Entropy of 3d BTZ black hole 
5.2.1 BTZ black hole geometry 

The black hole solution in three-dimensional gravity with negative cosmological constant was first 
obtained in [6] (see also for the global analysis of the solution). We start with the black hole 
metric written in a form that makes it similar to the four-dimensional Kerr metric. Since we are 
interested in its thermodynamic behaviour, we write the metric in the Euclidean form: 

ds^ = f{r)dT^ + r\r)dr^ + r^{d(j) + N{r)dTf , (224) 

where the metric functions /(r) and N[r) read 

Sir) 4 - ^; - .» ^ • - -v^ (^^^' 

and we use the notation 

-+-^(l + \A + (%)- KP-^Jl + (%-!) (226) 
2 V mi 2 s ml 

Obviously one has that r+jr-j = jl. The coordinate (j) in (|224p is assumed to be periodic with 
period 27r. 

In order to transform the metric (j224[) to Lorentzian singnature we need to make the analytic 
transformation t it, j — > —ij so that 

where and are the values in the Lorentzian space-time. These are the respective radii of the 
outer and inner horizons of the Lorentzian black hole in (2 -|- 1) dimensions. Therefore we must 
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always apply the transformation (|227|) after carrying out all calculations in the Euclidean geometry 
in order to obtain the result for the Lorentzian black hole. The Lorentzian version of the metric 
(j224l) describes a black hole with mass m and angular momentum J = 2j. The outer horizon is 
located at r — r+ , the respective inverse Hawking temperature is 

/'(r+) r+ + |r_|^ 

In the (r, r) sector of the metric (|224l) there is no conical singularity at the horizon if the Euclidean 

time r is periodic with period Ph- The horizon S is a one-dimensional space with metric ds^ = 

Pdtp"^ , where ij) = -j-4' ~ '^i^^''' is a natural coordinate on the horizon. 

The BTZ space is obtained from the three-dimensional maximally symmetric hyperbolic space 
(sometimes called the global Euclidean anti-de Sitter space) by making certain identifications. 

In order to see this one may use the coordinate transformation 

""h-V-C^^)'". (229) 

In new coordinates (p, 0^ ip) the BTZ metric takes the form 

ds^ = f {dp^ + cosh^ pdij^ + sinh^ pdO^) , (230) 

which is the metric on the hyperbolic space H^. In this metric the BTZ geometry is defined by 
identifications 
i). 0^e + 2TT 

a). 0^0 + 2Tr^-^, V^V^ + Stt^. 
The outer horizon r = in the coordinate system (p, 6, ip) is located at p ~ and ip is the 
angular coordinate on the horizon. Notice that the geodesic distance a between two points with 
coordinates {p,ip,9) and {p,tp' ,9') is 

• T,2 T,2 . , 2 -0 - "0' • , 2 . 2 6' - 61' 
smh — = cosh p smh — ^ h smh p sm — ^ — . (231) 

5.2.2 Heat kernel on regular BTZ geometry 

Consider a scalar field with the operator T) — — (V^ + C/^^)- The maximally symmetric constant 
curvature space is a nice example of a curved space in which the heat equation [ds+'D)K{x, x' , s) — 
has a simple, exact, solution. The heat kernel in this case is a function of the geodesic distance 
a between two points x and x' . On the global space Hj, one finds 

KHA<y.s)^ ^3,^ -wS/^ ^'^'"^' (232) 
{Airsyi'^ smh((T//) 

where /i = 1 — ^. The regular BTZ geometry is defined by identifications i) and ii) defined above. 
As is seen from (|231l) the geodesic distance and the heat kernel (|232l) . expressed in coordinates 
{p,il),9), are automatically invariant under indentification i). It remains thus to maintain the 
identification ii). This is done by summing over images 

Kbtz{x, x',s)= Kh, (p, p', V - 0' + 27r^n, 9 - 9' + 2^^n) . (233) 

n— — oo 

Using the path integral representation of the heat kernel we would say that the n = term in 
Ij233p is due to the direct way of connecting points x and x' in the path integral. On the other 
hand, the n ^ terms are due to uncontractible winding paths that go n times around the circle. 
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5.2.3 Heat kernel on conical BTZ geometry 

The conical BTZ geometry wtiicti is relevant to the entanglement entropy calculation is obtained 
from global hyperbolic space H3 by the replacing the identification i) as follows 

and not changing the identification ii). For a ^ 1 this Euclidean space has a conical singularity 
at the horizon [p — Q). The heat kernel on the conical BTZ geometry is constructed via the heat 
kernel (|233l) on the regular BTZ space by means of the Sommerfeld formula (1^^ 

Kbtz^{x,x',s) = Kbtz{x,x',s) + I cot ^ Kbtz{0 - 0' + w, s) dw , (234) 

47ra Jp Za 



where Kbtz is the heat kernel ([233]). The contour F is defined in (|22|) . 

For the trace of the heat kernel (|234p one finds |171j after computing by residues the contour 
integral 

1 tKbtz^ = — 7^ h —(2TTa)C2[a) s 



e-^^' A+ ^ sinh ^ e- 

+^''(4^ij^~r ' ^ sinh AV^„ (sinh^ M:. + sinh^ ^) ' 

where s = s/P, 7^ = A-n/l and A-0„ = A+n/l = 27rr+ and A^ — 27rr_). Notice that 
we have already made the analytical continuation to the values of and r_ in the Lorentzian 
geometry. 

5.2.4 The entropy 

When the trace of the heat kernel on the conical geometry is known one may compute the entan- 
glement entropy by using the replica trick. The entropy then is the sum of UV divergent and UV 
finite parts jl71J 

Sent = Sdiv + Sfin , (236) 

where the UV divergent part is 



This divergence is renormalized by the standard renormalization of Newton's constant 

^ 1 1 1 (238) 



in the three dimensional gravitational action. 
The UV finite part in the entropy is 



S fin — ^ ^ ^n J 



n=l 



(1 + A+n cothA+Ti 



2n (cosh^+n — cosh^_n) 
{A^n sinh A-^-n — A_n sinh A^n) 
(cosh A+n — cosh A^n) 



) , (239) 
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where A± — 2Trr±/l. 

After the renormahzation of Newton's constant the complete entropy of the BTZ black hole, 
Sbtz = Sbh + Senti IS a rather complicated function of the area of inner and outer horizons. 
Approximating in (I239P the infinite sum by an integral one finds [171] 

S = + / s{x) dx , 

1 g-v^^ / (2; sinh a: — /ex sinh fca;) \ f240) 

2x (cosh X — cosh kx) \ (cosh x — cosh kx) J ' 

where k = A^/A+. The second term in the right hand side of (|240p can be considered to be the 
one-loop quantum (UV-finitc) correction to the classical entropy of black hole. 

For large enough = -p >> 1 the integral in (|240l) goes to zero exponentially and we 
have the classical Bekenstein- Hawking formula for entropy. On the other hand, for small A_^, the 
integral in (|240p behaves logarithmically so that one has [171] 



_A+ JJIA+ 1/1+ 4+ 
Sbtz - + — — + «(( — ) ) • (241) 

This logarithmic behavior for small values of A+ (provided the ratio k = A^/A+ is fixed) is 
universal and independent of the constant ^ (or /i) in the field operator and the area of the inner 
horizon ( A_ ) of the black hole. Hence the rotation parameter J enters (j24ip only via the area A+ 
of the outer horizon. 

The other interesting feature of the entropy (j240p is that it always develops a minimum which 
is a solution to the equation 

I ,A+ 



AG '^^^ ■ 

This black hole of minimal entropy may be interesting in the context of the final stage of the 
Hawking evaporation in three dimensions. As follows from the analysis of Mann and Solodukhin 
[1711 . the minimum of the entropy occurs for a hole whose horizon area is of the Planck length, 
A+ ^ IpL (in threee dimensions IpL ^ Gren)- 



5.3 Entropy of d-dimensional extreme black holes 

The extremal black holes play a special role in gravitational theory. These black holes are character- 
ized by vanishing Hawking temperature Th which means that in the metric (jl65p the near-horizon 
expansion in the metric function g{r) starts with the quadratic term (r — r+)^. Topologically, the 
true extremal geometry is different from the non-extremal one. Near the horizon the non-extremal 
static geometry looks like a product of a two-dimensional disk (in the plane (r, r)) and a {d — 2)- 
dimensional sphere. The horizon then is the center in the polar coordinate system on the disk. 
Contrary to this, an extremal geometry in the near-horizon limit is a product of a two-dimensional 
cylinder and a (d — 2)-dimensional sphere. Thus, the horizon in the extremal case is just another 
boundary rather than a regular inner point as in the non-extremal geometry. However, one may 
consider a certain limiting procedure in which one approaches the extremal case staying all the 
time in the class of non-extremal geometries. This limiting procedure is what we shall call the 
"extremal limit" . A concrete procedure of this type was suggested by Zaslavsky [225] . One con- 
siders a sequence of non-extreme black holes in a cavity at r = and finds that there exists a 
set of data (r+, r^, r_) such that the limit ^ — 1, ^ — > 1 is well-defined. Even if one may have 
started with a rather general non-extremal metric the limiting geometry is characterized by very 
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few parameters. In this sense one may talk about "universality" of the extremal limit. In fact, in 
the most interesting (and tractable) case the limiting geometry is the product of two-dimensional 
hyperbolic space H2 with the (d — 2)-dimensional sphere. Since the limiting geometry belongs to 
the non-extreme class its classical entropy is proportional to the horizon area in accord with the 
Bekenstein-Hawking formula. The entanglement entropy of the limiting geometry then is a one- 
loop quantum correction to the classical result. The universality we have just mentioned suggests 
that this correction possesses a universal behavior in the extreme limit and, since the limiting 
geometry is rather simple, the limiting entropy can be found explicitly. The latter was indeed 
shown by Mann and Solodukhin in |172| . 



5.3.1 Universal extremal limit 

Consider a static spherically-symmetric metric in the following form 

ds^ = g{r)dT^ + ^dr^ + r^dujj_2 , (243) 

where dw^_2 ^^"^ metric on the {d — 2)-dimensional unit sphere, describing a non-extreme hole 
with an outer horizon located at r — r+. The analysis can be made for a more general metric, in 
which grr 7^ drr-: the limiting geometry however is the simplest in the case we consider in (j243[) . 
The function g{r) in (|243p can be expanded as follows 

g{r) = a{r - r+) + b{r - r+f + 0((r - r+f) . (244) 

It is convenient to consider the geodesic distance I — J g^^^^dr as a radial coordinate. Retaining 
the first two terms in (|244p . we find, for r > r+, that 

(r-r+) = -smh (^-) , 

g{lb'/') = {^)sinh\lb'/'). (245) 

In order to avoid the appearance of a conical singularity at r = r+, the Euclidean time r in p43|) 
must be compactified with period At: /a which goes to infinity in the extreme limit a — >■ 0. However, 
rescaling r — > = Ta/2 yields a new variable having period 2tt. Then, taking into account (|245l) 
one finds for the metric (|243p 

ds^ = i (sin^ xd(j)^ + dx^) + sinh^ i^fdujl-2 , (246) 

where we have introduced the variable x = Ib^/"^. To obtain the extremal limit one just takes 
a — >■ 0. The limiting geometry 

ds^ = i (sin^ xd<l)^ + dx^) + r\dujl_^ (247) 

is that of the direct product of a 2-dimensional space and a (d — 2)-sphere and is characterized by 
a pair of dimensional parameters and r+. The parameter r+ sets the radius of the (d — 2)- 

dimensional sphere while the parameter b~^/'^ is the curvature radius for the (a;, (j)) 2-space. Clearly, 
this two-dimensional space is the negative constant curvature space 7?2- This is the universality 
we mentioned above: although the non-extreme geometry is in general described by an infinite 
number of parameters associated with the determining function g{r) the geometry in the extreme 
limit depends only on two parameters b and r+. Note that the coordinate r is inadequate for 
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describing the extremal limit ()247p since the coordinate transformation (I245|) is singular when 
a — 0. The limiting metric (I247P is characterized by a finite temperature, determined by the 27r 
periodicity in angular coordinate 4>. 

The limiting geometry (|247l) is that of a direct product H2 x ^2 of 2d hyperbolic space H2 
with radius I — hr^l"^ and a 2d sphere S2 with radius = r+. It is worth noting that the limiting 
geometry (|247p precisely merges near the horizon with the geometry of the original metric ()243|) in 
the sense that all the curvature tensors for both metrics coincide. This is in contrast with, say, the 
situation in which the Rindler metric is considered to approximate the geometry of a non-extreme 
black hole: the curvatures of both spaces do not merge in general. 

For a special type of extremal black hole, when I = li, the limiting geometry is characterized 
by just one dimensionful parameter. This is the case for the Reissner-Nordstrom black hole in four 
dimensions. The limiting extreme geometry in this case is the well-known Bertotti-Robinson space 
characterized by just one parameter r_|_. This space has remarkable properties in the context of 
supergravity theory that are not the subject of the present review. 



5.3.2 Entanglement entropy in the extremal limit 

Consider now a scalar field propagating on the background of the limiting geometry (|247p and 
described by the operator 

I?=-(V2+X) , X = (248) 

where is the Ricci scalar. For the metric (|247p characterized by two dimensionful parameters 
I and h, one has that R(^d} — — 2/Z^ + {d — 2){d — For a d-dimensional conformally coupled 

4(d-l)- 



scalar field we have ^ = A^ri-^ \ ■ '^^^^ 



(d-2)(rf-4) (d-2) 
X,onf = +2(d^^^/^ (^^^^ 

The calculation of the respective entanglement entropy goes along the same lines as before. 
First, one allows the coordinate 0, which plays the role of the Euclidean time, to have period 27ra. 
For a 7^ 1 the metric (|247l) then describes the space E"' = H2 x 5^-2, where H2 is the hyperbolic 
space coinciding with H2 everywhere except the point x ^ where it has a conical singularity with 
an angular deficit S = 2tt{1 — a). The heat kernel of the Laplace operator on is given by 
the product 

Ke'-{z,z',s) = Kh»{x,x' ,s) Ks^_^{9,e' ,Lp,ip' ,s) 

where Kh^ and Ks^ are the heat kernels, of the Laplace operator on respectively H2 and Sd-2- 
The effective action reads 



1 d 

Weff[E'^] = ~-J^^ -TrKH^TrKs,_,e''' , (250) 

where e is a UV cut-off. On spaces with constant curvature the heat kernel function is known 
explicitly |33]. In particular, on a 2d space H2 of negative constant curvature the heat kernel has 
the following integral representation: 

1 \/2p"^/'* r°° dvve^y^/'^^ 
K^^(z,z',c) = ^_ / "yy*^ ^ (251) 

^ ' ' P (47rs)'^/^ Vcoshy — cosher ' 



where s — sl~'^. In equation (|251l) a is the geodesic distance between the points on H2- Between 
two points (a;, (/)) and (x, -I- A(/)) the geodesic distance is given by sinh^ ^ = sinh^oisin^ ^ . 
The heat kernel on the conical hyperbolic space H2 can be obtained from (I25ip by applying the 
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Sommerfeld formula (|22|) . Skipping the technical details, available in |172) . let us just quote the 
result for the trace 

— s/4 

TtKh^ = aTrifff, + (1 - q)-^---^:^^ (s) + 0(1 - a)^ , 

coshy 

' sinh^ y sinh 2y 



knis) = / ^2/^:^(1 - —TTtK' " ' (252) 



where s = s/P. 

Let us denote Qd-2{s) — TtKs^_2{s) the trace of the heat kernel of Laplace operator — on 
{d — 2)-dimensional sphere of unite radius. The entanglement entropy in the extremal limit then 
takes the form [172], [205] 



/ 4ifc/f(s)e,_2(4)e-^/V^'' . (253) 

J ^2/12 s-y^ If 



The function fc^f(s) has the following small-s expansion 

, 1 17 2 29 3 1181 4 1393481 . 763967 g ^ /or.^ 

knis) = y/irsi s H s"* H s'^ H s*" + ..) (254) 

-"^ ^ ^3 20 1120 4480 337920 615014400 447283200 ' 

The trace of the heat kernel on a sphere is known in some implicit form. For our purposes however 
a representation in a form of an expansion is more useful, 

Q / OO \ 

e'^-2(.) = J^^^^hm I 1 + (rf - 2)(d ~ 3) ^ a.„." ) , (255) 



n=l 



where = r^(|rf_i)/2) ^'^^ area of a unit radius sphere Sd-i- The first few coefficients in this 

expansion can be calculated using the results collected in [218] , 

1 (5d2-27d + 40) (35^* - 392^3 + 1699^2 - 3322d + 2520) 
= 6 ' ^ 360 ' ^ i5360 ' ^'^^^ 

We shall consider some particular cases. 
d=4. The entanglement entropy in the extreme limit is 



'5.=4 = J^ + solny + s(^), .o = ^(6e-l) + ^|(l-50, (257) 



where s(y) is UV finite part of the entropy. For minimal coupling (^ = 0) this result was obtained 
in |172j . The first term in (I257P is proportional to the horizon area A = 47r^^ while the second 
term is a logarithmic correction to the area law. For conformal coupling ^ = 1/6 the logarithmic 
term is 
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d=5. The entropy is 



^.=5 + ~ 5P) + lOm' - /?))^ + -(y) . (259) 

To simplify the expressions in higher dimensions we consider only the case of the conformal coupling 

c _ d~2 
? - 4(d-l)- 

The entropy takes the form: 
d=6. 



d=7. 



d=8. 



y/^q Ty/^ (25P - 32/^)^3 
384e5 ^ 34560 

(70592/f - 109760/^/2 + 40635/-^)/! h 



1935360 l^e ' / 



- + Kt)- (261) 



S,=s^'4+'-^+'4 + '^0 In e + ,s(^) , (262) 
e° e"' i 

_ _ 1 /f(-209/2 + 160/2) 

~ 2160 ' ~ 75600 /6 ' 

_ 1 /^(8753/t - 13376/2/2 + 4875/^) 

~ 352800 P ' 

1 (1102263/f- 2520864/2/4 + 1833975/^/? - 413120/6) 

- 11113200 ■ 



Two examples of the extreme geometry are of particular interest. 



Entanglement entropy of the round sphere in Minkowski spacetime. Consider a sphere 
of radius R in flat Minkowski spacetime. One can choose a spherical coordinate system (r, r, 9^) so 
that the surface S is defined as r = and r = R, and variables 6*', i = 1, .., d — 2 are the angular 
coordinates on E. The d- metric reads 

ds2 = dr^ + dr^ + r^dujj_2 , (264) 

where duj'^_2 is metric on {d — 2) sphere of unite radius. Metric (|264p is conformal to the metric 

dsLt - ^(dr2 + rfr2) + R'dujj_, , (265) 

which describes the product of two-dimensional hyperbolic space H2 with coordinates (r, r) and 
the sphere 5'd-2- Note that both spaces, H2 and Sd~2, have the same radius R. Metric (|265l) 
describes the spacetime which appears in the extremal limit of a d-dimensional static black hole. 
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So 



90 



2881 
756 



1569275563 
1111320 



Table 1: Coefficients of the logaritfimic term in the entanglement entropy of extreme Reissner- 
Nordstrom black hole. 



In the hyperbolic space H2 we can choose a polar coordinate system (p, (p) with its center at point 

r = R, 



R 



R sinh p sin ( 



cosh p — sinh p cos i 

(for small p one has that r = R + p cos , 
that the metric takes the form 



(266) 



i cosh p — sinh p cos (p 

T — p sin (j) as in the polar system in flat spacetime) so 



dsl^^ = R^{dp^ + sinh^ 



(267) 



In this coordinate system the surface S is defined by the condition p — 0. In the entanglement 
entropy of a conformally coupled scalar field the logarithmic term sq is conformal invariant. There- 
fore, it is the same |205j for the entropy of a round sphere of radius R in Minkowski spacetime 
and in the extreme limiting geometry (j267|) . In various dimensions sq can be obtained from the 
results (|258|) - ()262p by setting I = li = R. One finds So = ^ in d = 4, sq = "yge in = 6 and 
So = -^j^gi^QQ in d = 8. For c? > 4 the logarithmic term in the entropy of a round sphere has been 
calculated by Casini and Huerta [38] directly in Minkowski spacetime. They have obtained sq in 
all even dimensions up to d = 14. Subsequently, Dowker [71] has extended this result to d = 16 
and d = 18. In arbitrary dimension d the logarithmic term so can be expressed in terms of Bernulli 
numbers as is shown in 1381 and 1711 



Entanglement entropy of the extreme Reissner-Nordstrom black hole in d dimensions. 

As was shown by Myers and Perry [179] a generalization of the Reissner-Nordstrom solution to 
higher dimension d > 4 is given by (|243p with 



,d-3 



,d-3^^^d-3 



j,2(ci-3) 



(268) 



In the extreme limit 



Expanding (|268p near the horizon one finds, in this limit, that 



6 = (d — 3)^/r^ . Thus this extreme geometry is characterized by values of radii I = r_|_/(d — 3) and 
li = r^. In dimension d = 4 wc have ^ = Zi as in the case considered above. In dimension d > 4 
the two radii are different, I ^ li. For the conformal coupling the values of logarithmic term sg in 
various dimensions are presented in Table 1. 



6 Logarithmic term in the entropy of generic conformal 
field theory 

As we have already seen, in even dimensions, there typically appears a logarithmic term in the 
entanglement entropy. This term is universal in the sense that it does not depend on the scheme 
which is used to regularize the UV divergences. In conformal field theories the logarithmic terms 
in the entropy are closely related to the conformal anomaly. In this section we discuss in detail this 
aspect and formulate precisely the relation between entanglement entropy and conformal anomalies. 
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Consider a conformal field theory in d spacetime dimensions. As we have discussed throughout 
this review the most efficient way to calculate the corresponding entanglement entropy for a non- 
extremal black hole is to introduce a small angle deficit S = 27r(l — a) at the horizon surface S, 
compute the effective action W{E°') on a manifold E" with a singular surface and then apply 
the replica formula S — {ada — i)\a=iW{E") and obtain from it the entanglement entropy. In d 
dimensions the effective action has the general structure 

WcME-)^^-^-..-^-..-a,,,lne + w{gi-^) . (269) 

The logarithmic term in this expansion appears only if dimension d is even. Thus only even d 
will be considered in this section. The terms a^;, ad-2j-- representing the power UV divergences, 
are not universal, while the term 0^/2 is universal and is determined by the integrated conformal 
anomaly. The term w{g'-°'')) is the UV finite part of the effective action. Under a global rescaling 
of the metric on E" , g^"^ — \'^g'^"\ one has 

u;(AV"^) =w(5^"^) + ad/2lnA . (270) 

An important property of the expansion (|269p for a quantum field theory which classically 
is conformally invariant is that the logarithmic term aj;/2 is conformal invariant (see [218] and 
references therein), 

ad/2[e^^"g] = ad/2 [5] • 

On a manifold with a conical singularity at the surface E, the coefficients ad-2n have a bulk part 
and a surface part. To first order in (1 — a) one finds that 

«d-2„(i?") - aa^d'^LiE) + (1 - a)a^-2n + 0{1 - (271) 

For n — d/2 the oefficients a^^^{E) and 0^/2 (S) are respectively the integrated bulk and surface 
conformal anomalies. The bulk and surface terms are independently invariant under conformal 
transformation g — >■ e~'^^ g, 

«d/2'(e-'"5) - 4f2{9) and a^„2{e-'-g) = a%M . (272) 
Applying the replica formula one obtains the entanglement entropy 

So = af/2 : Sd-2n = , n = 1, .., d/2 - 1 

s{\^g)^s{g)-a^,^\n\ . (273) 

For a regular manifold each term ad-2n is the integral of a polynomial of degree n in the Riemann 
curvature. Respectively, the surface terms af_2„ is the integral over the singular surface S of poly- 
nomial of degree n — 1 in the Riemann curvature and its projections onto the subspace orthogonal 
to surface S. The concrete structure of the polynomials depends on the dimension d. 

6.1 Logarithmic terms in 4-dimensional conformal field theory 

In four dimensions the bulk conformal anomaly is a combination of two terms, the topological 
Euler term and the square of the Wcyl tensor, 

£^(4) = ^ / {RaPfiiyR"'^^'' — '^Rfj.uR^" + R^) , 
04 J E 

1(4) = ^(i?a/3M.i?"^^'' - 2i?^.i?^'^ + ii?^) . (274) 
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These are, respectively, the conformal anomalies of type A and B. In a theory with Ug particles of 
spin s one finds [75] (the contributions of fields of spin 3/2 and 2 can be obtained from table 2 on 
p.l80 of the book of BirreU and Davies P]) 

A = -^{na + llni/2 + 62?ii + 0713/2 + 0^2) , 

1 233 424 

B = 1^77^ ("-0 + 6ni/2 + 12rii 7-^3/2 + ^"2) • (275) 

307r^ ' 6 ' 3 

The surface contribution to the conformal anomaly can be calculated directly by, for example, the 
heat kernel method as in |103j . The direct computation although straightforward is technically 
involved. One has however a short cut: there is a precise balance, observed in [197 and 
between the bulk and surface anomalies, this balance is such that, to first order in (I — a), one can 
take a2{E") = 03 ""^(-E") + 0(1 — a)^ and use for the Riemann tensor of E" the representation as a 
sum of regular and singular (proportional to a delta-function concentrated on surface E) parts. The 
precise expressions are given in [111) . [112) . This representation, however, is obtained under the 
assumption that the surface S is a stationary point of an Abelian isometry and thus has vanishing 
extrinsic curvature. Under this assumption one finds that |112j (see also [187] ) 

a2iE") = aaY^^{E) + (1 - a)a| + 0(1 - af , 



= Aa^ + Ba% , 

— J [Rijij — 2R^i + R) 



IT, 

{Rijij — Rii + z^R) , (276) 



s 



3 



where Rijij = j?affp,.n°n^rtf nj, Raa = Rapn'^n'^. 

Each surface term in (|276ji is invariant under a sub-class of conformal transformations, g — > 
e~^"(7, such that the normal derivatives of uj vanish on surface E. The surface term due to the 
bulk Euler number is, moreover, a topological invariant: using the Gauss-Codazzi equation 

R^Rt. + 2Ru - R^j^J -k'k' + Trfc2 , (277) 

where i?s is the intrinsic Ricci scalar of the surface and fc^^ is the extrinsic curvature, and in the 
assumption of vanishing extrinsic curvature the term, as shown in [112J . is proportional to the 
Euler number of the 2d surface S, 

al = ^j^RT , (278) 

where i?s is intrinsic curvature of S. 

For completeness we note that this result can be generalized to an arbitrary codimension 2 
surface in 4-dimensional spacetime. The conformal transformation then is generalized to any 
function uj with non- vanishing normal derivative at E. The terms with the normal derivatives of w 
in the conformal transformation of then can be cancelled by adding the quadratic combinations 
of extrinsic curvature, Tr/c^ and kaka- The analysis presented by Solodukhin [203] (this analysis is 
based on an earlier consideration by Dowker |75] ) results in the following expressions 

= Aa^ + Sa| , 

aA = l [ (Rm - + R- TrA:^ + k,h) = I [ Rs , 

TT f 1 1 

a| = - J JJR^3^J - R^^ + 3^ " (Trfc' - l^hh)) . (279) 



59 



This is the most general form of the fogarithmic term in the entanglement entropy in four spacetime 
dimensions. 

Thus, as follows from equation (I274p . the logarithmic term in the entanglement entropy of black 
hole in four dimensiosn is 

4-^=4) = A^JjR,j,, - 2Ru + R)~B^ JjR,,,j - Ru + \r) . (280) 

For conformal fields of various spin the values of A and B are presented in (|275p . 
Consider some particular examples. 

Extreme static geometry. For an extreme geometry which has the structure of the product 
H2 X S2 and characterized by two dimensionful parameters I (radius of H2) and li (radius of S2) 
the logarithmic term in the entropy 

sr = An'-^{l-f^) (281) 

is determined by both the anomalies of type A and B. In the case of the extreme Reissner-Nordstrom 
black hole one has I — li and the logarithmic term (I28ip is determined only by the anomaly of 
type A. For a conformal scalar field one has that A = B /3 = I/QOtt^ and this equation reduces to 
(|258p . As we already discussed, the geometry H2 x S2 for I = li is conformal to flat 4-dimensional 
space. Thus the Weyl tensor vanishes in this case as does its projection to the subspace orthogonal 
to horizon 82- That is why the type B anomaly does not contribute in this case to the logarithmic 
term. 

The Schwarzschild black hole. In this case the background is Ricci flat and the logarithmic 
term is determined by the difference of A and B, 

sf' ^{A- B)Tr^ . (282) 

The same is true for any Ricci flat metric. For a conformal scalar field equation (I282p reduces 
to (|55)) . For a scalar field the relation of the logarithmic term in the entropy and the conformal 
anomaly was discussed by Fursaev |105J . The logarithmic term vanishes ii A = B. In this case the 
Riemann tensor does not appear in the conformal anomaly (|274p so that the anomaly vanishes if 
the metric is Ricci fiat. In particular, the relation A = B can be found from the Af = A super- 
conformal gauge theory, dual to supergravity on AdS's, according to the AdS/CFT correspondence 
of Maldacena |167) . The conformal anomaly in this theory was calculated in [131) . 

Non-extreme and extreme Kerr black hole. For the Kerr black hole (characterized by mass 
m and rotation a) the logarithmic term does not depend on the parameter of the rotation and it 
takes the same form 

^Kerr ^ _ B)7r^ (283) 

as in the case of the Schwarschild metric. In the extreme limit a = m the logarithmic term takes 
same value (|283p . 
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6.2 Logarithmic terms in 6-dimensional conformal field theory 

In six dimensions, omitting the total derivative terms, the conformal anomaly is a combination of 
four different conformal invariants [T7] 

abuik^f (^B^i^ + B2I2 + B3I3 + AEe) , (284) 
where Ii , I2 and are cubic in the Weyl tensor 



J /^imnj kl j kl mn ij 

^1 — ^kijl'-^ L-m „ , -(2 — ^ij ^kl ^mn i 

6 

5 

and is the Euler density (|50)) 



^3 = {V^5] + 4i?; - . (285) 



30727r3 i'ii>2---^ I'si'e • l,^ou; 

As was shown in [17_ in a free conformal field theory with ng scalars, 711/2 Dirac fermions and 
2-form fields one has thalF°l 

, 8-3!, 5 191 221 , 

A = ^{-—nQ-—n^/2~—nB), (287) 
1/28 896 8008 , 



1,5 ^„ 2378 



^2 = 773T^(-:7"o + 32ni/2 + ^^nB), (289) 



-^3 = 7r43^(-2"o-40ni/2- 180ns). (290) 
(47r)'^7! 

Applying the formulas ((55|) to /i, /2 and /a and using the relation (|62|) for the Euler number 
one finds for the logarithmic term in the entanglement entropy of 4-dimensional surface S in a 
6-dimensional conformal field theory 

S't^ = Bi SBl + B2 SB2 + S3 SB3 + ^ SA , (291) 

where 

= XP] , (292) 

is the Euler number of the surface S, and 

SBl = 6^(C^"^'CJ^' p - C^^P^C^' p - \G"^^>'C\f,^ + Ic^^^'-C^;,^,) , (293) 

,ss2 = 6^(2C'^"^C„;-'' - C'^-^^C^f,^ + ic"''^''a/3^,) , (294) 

SB3 = 87r(V2C*J''J + 4i?*„C"^'*^' - i?„^C"*''^ - ^i?C*^''J' + C"„^^C""^'' - ^C"''''''Ca;3M^) , (295) 

where tensors with Latin indices are obtained by contraction with components of normal vectors 
nj,, i = 1,2. Note that in (|295|) we used for brevity the notation V^C^-?*^ = n'^nj^n'-^nlV'^C"'^'"' . 
Equations p93)) . p94)) . p95|) agree with result obtained in fT34l. 

Let's consider some examples. 



Note that the coefficient be of paper | 17| is related to 03 as feg = —13. 
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6-dimensional Schwarzschild black hole. The 6-dimensional generalization of the Schwarzschild 
solution is [179] 



r2 



ds^ = g{T)dT^ + g-\r)dr^ + r^duj^^ , g{r) = 1 - ^ , (296) 

where is metric of unit 4-sphere. The area of horizon is yl+ = ^^r^^. The Euler number of 
the horizon x[S'4] = 2. This metric is Ricci flat so that only the Riemann tensor contributes to the 
Weyl tensor. The logarithmic term in this case is 



Sch 



167r^(-51Bi + 156^2 - I92S3) + 2A . (297) 



It is interesting to note that this term vanishes in the case of the interacting (2, 0) conformal theory 
which is dual to supergravity on AdSj. Indeed in this case one has [131) . [T7] 

35 

a = -y , 61 = 1680 , 62 = 420 , 63 = -140 (298) 

so that s[,^'°^ = 0. This is as expected. The Riemann tensor does not appear in the conformal 
anomaly of the strongly interacting (2, 0) theory so that the anomaly identically vanishes if the 
spacetime is Ricci flat. This property is not valid in the case of the free (2, 0) tensor multiplet [17] 
so that the logarithmic term of the free multiplet is non- vanishing. 



Conformally flat extreme geometry. In conformally flat spacetime the Weyl tensor Cap^iv — 
so that terms (|293p . (I294p and (|295l) identically vanish. The logarithmic term (|291l) then is 
determined by the anomaly of type A only. In particular this is the case for the extreme geometry 
H2 X 5*4 with equal radii Z = Zi of two components. One has 

sf;^* = 2A (299) 

for this extreme geometry. This geometry is conformal to flat spacetime and the logarithmic 
term ([299]) is the same as for the entanglement entropy of a round sphere in flat 6-dimensional 
spacetime. This generalizes the result discussed in section [5.3.21 for the entropy of a round sphere 
due to conformal scalar field. The result (I299p . as is shown in 09], |180j . generalizes to a spherical 
entangling surface in a conformally fiat spacetime of any even dimension. 



6.3 Why might logarithmic terms in the entropy be interesting? 

By a logarithmic term we mean both the logarithmically UV divergent term in the entropy and 
the UV finite term which scales logarithmically with respect to the characteristic size of the black 
hole. As we have seen, these terms are identical. However, after the UV divergences in the entropy 
have been renormalized it is the UV finite term, which scales logarithmically, that will interest us 
here. 

i) First of all, the logarithmic terms are universal and do not depend on the way the entropy was 
calculated and on the scheme in which the UV divergences are regularized. This is in contrast with 
the power UV divergences in the entropy that depend both on the calculation procedure and on 
the regularization scheme. 

it) The logarithmic terms are related to the conformal anomalies. As the conformal anomalies play 
an important role in the modern theoretical models, any new manifestation of the anomalies merits 
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of our special attention. This may be even more important in view of ideas that the conformal 
symmetry may play a more fundamental role in Nature than is usually thought. As is advocated 
by 't Hooft in a number of recent papers |217) . [216) . |215j a crucial ingredient for understanding 
Hawking radiation and entropy is to realize that gravity itself is a spontaneously broken conformal 
theory. 

iii) For a large class of extremal black hole solutions, which arise in supergravity theories considered 
as low energy approximation of string theory, there exists a microscopic calculation of the entropy. 
This calculation requires a certain amount of unbroken supersymmetry, so that the black holes in 
question are the BPS solutions and uses the conformal field theory tools, such as the Cardy formula. 
The Cardy formula predicts certain logarithmic corrections to the entropy (these corrections are 
discussed, in particular, in |35j and |200) ). One may worry whether exactly same corrections are 
reproduced in the macroscopic, field theoretical, computation of the entropy. This aspect was 
studied recently in [8] for black holes in A/" = 4 supergravity and at least some partial (for the 
entropy due to matter multiplet of the supergravity) agreement with the microscopic calculation 
has been indeed observed. 

iv) Speaking about the already renormalized entropy of black holes and taking into account the 
backreaction of the quantum matter on the geometry, the black hole entropy can be represented as 
a series expansion in powers of Newton's constant, G/r^ (the quantity l/?'^, where Vg — 2GM is 
the size of black hole, is the scale of the curvature at the horizon; thus the ratio G/r^ measures the 
strength of gravitational self- interaction at the horizon) or, equivalently, in powers of M|,^/M^. 
In particular for the Schwarzschild black hole of mass M in four spacetime dimensions one finds 



S^Anj^ + a\nM + 0{^). (300) 



The logarithmic term is the only correction to the classical Bekenstein-Hawking entropy that is 
growing with the mass M. 

v) Although the logarithmic term is still negligibly small compared to the classical entropy for 
macroscopic black holes it becomes important for small black holes especially at the latest stage 
of the black hole evaporation. In particular, it manifests itself in a modification of the Hawking 

temperature as a function of mass M [105] . Indeed, neglecting the terms O(^^^r^) in the entropy 
([SUni one finds 

1/Th ^gj^^ SttM/M^^ + aAf-i (301) 

so that the Hawking temperature Th ~ M for small black holes. Depending on the sign of the 
coefficient a in (|300l) there can be two different scenarios. If cr < then the entropy S{M) as 
function of mass develops a minimum at some value of Mmin ~ Mp^. For this value of the mass 
the temperature p01[) becomes infinite. This is the final point (at least in this approximation) of 
the evaporation for black holes of mass M > Mmin- It is reached in finite time. Not worrying 
about exact numerical factors one has 

^ ^ -Tf,A+ ^ -Tt,M' (302) 

for the evaporation rate. For large black holes (Mq ^ Mmin) the evaporation time is Ibh ~ 
M^/Mp^. This evaporation time can be obtained by solving p02|) with the classical expression 
for the Hawking temperature Th ~ Mpj^/M, i.e. without the correction as in poip . Thus, if 
there is no logarithmic term in pOOp any black hole evaporates in finite time. If the correction 
term is present, it becomes important for Mg ~ Mmin- Then, assuming that Mmin ~ MpL, one 
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finds tsH ~ {Mq — Mpl)^ / Mp^ for the evaporation time. On the other hand, a black hole of 
mass Mq < Mmim if it exists, evaporates down to zero mass in infinite time. Similar behavior is 
valid for black holes of arbitrary mass Mq if ct > in equations pOO[) . poip . The evaporation 
rate considerably slows down for small black holes since the Hawking temperature Th M for 
small M. The black hole then evaporates to zero mass in infinite time. Asymptotically, for large 
time t, the mass of black hole decreases as M{t) ~ t~^/^ . Thus, the sub-Planckian black holes if 
(a < 0) and any black holes if (cr > 0) live much longer than one would have expected if one used 
the classical expression for the entropy and for the Hawking temperature. 

7 Holographic description of entanglement entropy of black 
hole 

A popular trend in modern fundamental physics is to reconsider various, sometimes very well 
known, phenomena from the point of view of holography. Holography is a rather general statement 
that the physics inside a spatial region can be understood by looking at a certain theory defined 
on the boundary of the region. This so-called holographic principle was first formulated by 't 
Hooft [214] and later generalized by Susskind |211'. A review of the holographic principle was 
given in [21]. A concrete realization of holography is the AdS/CFT correspondence |167| . [222] . 
[124] . According to this correspondence the theory of supegravity (more precisely string theory 
the low energy regime of which is described by a supegravity) in a (d + l)-dimensional anti-de 
Sitter spacetime (AdS) is equivalent to a quantum conformal field theory (CFT) defined on a 
d-dimensional boundary of the anti-de Sitter. There is a precise dictionary of how phenomena 
on one-side of the correspondence can be translated into phenomena on the other side. The 
correspondence has proved to be extremely useful, both for better understanding the gravitational 
physics and the quantum field theory. If c? = 4 then the CFT on the boundary is known to be a 
A/" = 4 superconformal gauge theory. This theory is strongly coupled and in many aspects resembles 
the QCD. Thus utilizing the correspondence one, in particular, may gain some information on how 
the theories of this type behave (for a review on the correspondence and its applications see [T]). 

One of the aspects of the AdS/CFT correspondence is geometrical. The boundary theory 
provides certain boundary conditions for the gravitational theory in the bulk so that one may 
decode the hologram: reconstruct the bulk spacetime from the boundary data. As was analyzed in 
|60j for this reconstruction the boundary data one has to specify consist on the boundary metric 
and the vacuum expectation of the stress-energy tensor of the boundary CFT. The details are 
presented in [B^, see also review |192| . 

Entanglement entropy is one of the fundamental quantities which characterize the boundary 
theory. One would think that it should have an interpretation within the AdS/CFT correspon- 
dence. This interpretation was suggested in 2006 by Ryu and Takayanagi [188] . [187] (for a review 
on this proposal see [184] ). This proposal is very interesting since it allows one to compute the 
entanglement entropy in a purely geometrical way (see also [109] ). 

7.1 Holographic proposal for entanglement entropy 

Let M be a (c? -I- l)-dimensional asymptotically anti-de Sitter spacetime. Its conformal boundary 
is a d-dimensional spacetime dM . On a slice of constant time t in dM one picks a closed {d — 2)- 
dimensional surface E and defines the entropy of entanglement with respect to S. Now, on the 
constant t slice of the (d -f l)-dimensional anti-de Sitter spacetime consider a (d — l)-dimensional 
minimal surface 7 such that its boundary in dM is the surface S, — S. According to the 
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proposal of |188j . [187J the following quantity 



(303) 



where ^(7) is the area of minimal surface 7 and Gd+i is Newton's constant in the ((i+l)-dimensional 
gravitational theory, is equal to the entanglement entropy one has calculated in the boundary 
conformal field theory. This holographic proposal has passed many tests and never failed. It has 
reproduced correctly the entropy in all those cases when it is known explicitly. In particular, in 
two spacetime dimensions (d = 2) it correctly reproduced pHt and pT|) for the entropy at finite 
size and at finite entropy respectively. In higher dimensions (d > 2) the area of minimal surface 
7 diverges when it is extended till dM. This is an important feature, typical of the AdS/CFT 
correspondence. In fact, instead of the conformal boundary dM one has to consider a regularized 
boundary dM^ located at a small distance e (measured in terms of some radial coordinate p). In 
the AdS / CFT correspondence the divergence in e has the interpretation of a UV divergence in the 
boundary quantum field theory. Considering the regularized surface 7 which extends to dM^, one 
finds that its area, to leading order in e, behaves as ^(7) ~ A(E)/e'^~^. Taking this behavior of 
the area one sees that the proposal pOSp correctly reproduces the proportionality of the entropy 
to the area of surface E and its dependence on the UV cutoff e. 

We note however that in certain situations the choice of the minimal surface E may not be 
unique. In particular, if the dividing surface E has several components or if the quantum field 
resides inside a cavity instead of being defined on an infinite space, there is more than one natural 
choice of the minimal surface 7. Different choices may correspond to different phases in the 
boundary theory |157j . 



7.2 Proposals for holographic entanglement entropy of black hole 

If one wants to generalize the proposal of Ryu and Takayanagi to black holes, the first step would 
be to find a (d + l)-dimensional asymptotically AdS metric which solves the Einstein equations 
in the bulk and whose conformal boundary describes a c?-dimensional black hole. To find such a 
metric explicitly may be a difficult task although some exact solutions are known. First of all, it 
is easy to construct an AdS metric which gives a de Sitter spacetime on the boundary. A de Sitter 
horizon is in many aspects similar to a black hole horizon. Entanglement entropy associated to a 
de Sitter horizon |128| . [139] has the same properties as the entropy of any other Killing horizon. 
In four dimensions {d = 3) an exact solution to the Einstein equations has been found in '84' that 
describes a 3c?-black hole on the boundary. In three dimensions {d = 2) an exact solution which 
describes a generic static two-dimensional black hole on the boundary has been found in [193] . On 
the other hand, the results of 60^ show that for any chosen metric on the boundary dM one can 
find, at least in a small region close to the boundary, an exact solution to Einstein equations with 
negative cosmological constant. Exact formulas are given in [60 . Thus, at least principally, it is 
not a problem to find an asymptotically AdS metric which describes a black hole on the boundary. 

The next question is how to choose the minimal surface 7. A proposal of Emparan |83j consists 
of the following. Suppose the metric on the boundary of asymptotically AdS spacetime describes 
a black hole with a Killing horizon at surface E. Presumably, the horizon on the boundary dM is 
extended to the bulk. The bulk horizon is characterized by vanishing extrinsic curvature and is a 
minimal {d — l)-dimensional surface. Thus, one can choose the bulk horizon to be that minimal 
surface 7, the area of which should appear in the holographic formula p03|) . In this construction 
the Killing horizon E is the only boundary of the minimal surface 7. This prescription is perfectly 
eligible if one computes the entanglement entropy of a black hole in infinite spacetime. In [53] it 
was applied to entropy of a black hole residing on the 2-brane in the 4d solution of ref.[Sl]. In 
ref.[139j a similar prescription is used to compute entanglement entropy of de Sitter horizon. 
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On the other hand, in certain situations it is interesting to consider a black hole residing inside 
a cavity, the so-called "black hole in a box". Then, as we have learned in two-dimensional case, 
the entanglement entropy will depend on the size L of the box so that, in the limit of large L, the 
entropy will have a thermal contribution proportional to volume of the box. This contribution is 
additional to the purely entanglement part which is due to presence of the horizon S. In order to 
reproduce this dependence on the size of the box one should use a different proposal. A relevant 
proposal was suggested by Solodukhin in [20 IJ . 

Let a d-dimensional spherically symmetric static black hole with horizon E lie on the regularized 
boundary (with regularization parameter e) of asymptotically anti-de Sitter space-time AdSd+i 
inside a spherical cavity of radius L. Consider a minimal d-surface T whose boundary is the 
union of S and E^. T has saddle points where the radial AdS coordinate has an extremum. 
By spherical symmetry the saddle points form a (d-2)-surface C with the geometry of a sphere. 
Consider the subset Vh of F whose boundary is the union of E and C. According to prescription 
of |201] . the quantity 

S = (304) 

is equal to the entanglement entropy of the black hole in the boundary of AdS. In particular, it 
gives the expected dependence of the entropy on the UV regulator e. The minimal surface 
"knows" about the existence of the other boundary E^. That is why (|304p reproduces correctly 
the dependence of the entropy on the size of the "box" . In |201) this proposal has been verified for 
d=2 and d = A. 

It should be noted that what as far as the UV divergent part of the entanglement entropy is 
concerned, the two proposals [53], [139] and [201] give the same result. This is due to the fact that 
the UV divergences come from that part of the minimal surface which approaches the boundary 
dM of the AdS spacetime. In both proposals this part of the surfaces 7 and F;i is the same. The 
difference thus is in the finite terms which are due to global properties of the minimal surface. 

From geometrical point of view the holographic calculation of the logarithmic term in the 
entanglement entropy is related to the surface anomalies studied by Graham and Witten [123] 
(this point is discussed in [190] ). 

7.3 Holographic entanglement entropy of 2d black holes 

In order to check the proposal p04p one needs a solution to the bulk Einstein equations that 
describes a black hole on the boundary of AdS. In three dimensions a solution of this type is 
known explicitly [193] . 



,2 dp^ 
ds = + - 
Ap^ p 




, (305) 



where the AdS radius is set to unity. At asymptotic infinity {p = 0) of the metric pOSp one has 
the 2d black hole metric 

dsla = f{x)dT^ + f-\x)dx^ , (306) 

where /(x) has simple zero in a: = x+. The cavity E^, is placed at a; = L so that < x < L. The 
regularity of the metric p06l) at the horizon a; = x+ requires that < r = 47r//'(a;+). 

Note, that pOSp is a vacuum solution of the Einstein equations for any function f{x). The regularity 
of the 3d metric p05l) requires that the constant h should be related to the Hawking temperature 
of the two-dimensional horizon by 6 = f'{x+). The geodesic F lies in the hypersurface of constant 
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time r. The induced metric on the hypersurface (p, x) has a constant curvature equal —2 for any 
function f{x) and is, thus, related by a coordinate transformation to the metric 



dsi 



(307) 



The exact relation between the two coordinate systems is 



1 



w 



z{x) 



/ l6/(x)-(6^-/^)p 
^16/(x) + (6-/i)V 



(16/(a;) + (fe-/i)2p) 



(308) 



where z(a;) = | ^^i)' '^'^^ equation for the geodesic in metric p07p is r = ^ — (w — . The 
geodesic length between two points lying on the geodesic with radial coordinates ri and r2 is 



7(M) = - 



In 



1 - Vl - C^r 



1 + Vl - C^r 



The saddle point of the geodesic is at — The constant C is determined from the condition 

that geodesic F joins points {x = x+,p = e^) and {x = L,p = e^) lying on the regularized (with 
regularization parameter e) boundary. In the limit of small e one finds that 



exp(— 6 



dx 

m 



The length of the geodesic joining point r-|_ corresponding to {x — x^,p = e) and the saddle 
point is thus 7(r;i) 



■In 



2 4 ■ 

Now, one has to take into account that, in the AdS/CFT correspondence, the value of Newton's 
constant in the bulk is related to the number of quantum fields living in the boundary DM. In the 
AdS^ / C FT2 case one has that ^r— — ^, where c is the central charge of boundary CFT. One thus 
finds the holographic entropy (|304l) 
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(309) 



dx 

W) 



{b-r) + lnf{L)-lnb^ 



where b = /'(x+), indeed coincides with the expression (j219p for the holographic entanglement 
entropy of the 2d black hole in conformal field theory. In particular, for large values of L the 
holographic formula for the entropy correctly reproduces the entropy of thermal gas Sth = ^ThL 
(we remind that Th = f {x+) / Att). This is a consequence of the choice of the minimal surface F/j 
made in the proposal p04|) . 



7.4 Holographic entanglement entropy of higher dimensional black holes 

In higher dimensions there is no known exact solution similar to p05|) . However, a solution in 
the form of a series expansion in p is available. In the rest of this section we focus on the case of 
boundary dimension 4. Then one finds [131] 

4p2 p 

]{x,p) = g(o)(a;) +ff(2)P + 5(4)P^ + /i(4)P^ Inp^ 



gij{x, p)dx^dx-' 



(310) 
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where g{Q)ij{x) is the boundary metric, coefficient |131j 

5(2)y = --^^Rtj - Q^gioyij) (311) 

is the iocai covariant function of boundary metric. Coefficient (7(4) is not expressed as a focal 
function of the boundary metric and is related to the stress-energy tensor of the boundary CFT 
|60] . which has an essentially nonlocal nature. /i(4) is a local, covariant, function of the boundary 
metric and is obtained as a variation of the integrated conformal anomaly with respect to the 
boundary metric [BD]- Its explicit form was computed in [BD]- ^(4) is a traceless tensor and in 
mathematics literature it is known as the "obstruction" tensor |122| . The explicit form of /i(4) or 
5(4) is not important if one wants to compute the UV divergence terms in the entropy. 

One may choose the boundary metric describing a static spherically symmetric black hole to 
take the form 

ds^ = f{r)dT^ + r\r)dr^ + r^{d0^ + sin^ Odcp^) (312) 

The minimal surface F lies in the hypersurface of the constant r of 5-dimensional space-time piO|) . 
The induced metric on the hypersurface takes the form 



,2 dp' 1 
4p2 p 



2 



F^^+R\de' + siu' 
fin 



(313) 



where functions F{r,p) — g^Q-^grr and R^{r,p) — ggg have p-expansion due to piO[) . The minimal 
surface F can be parameterized by (p, 9, 0). Instead of the radial coordinate r it is convenient to 
introduce the coordinate y = J dr/y/J so that, near the horizon, one has r = r+ + by'' /A + 0{y'^). 
The coordinate y measures the invariant distance along the radial direction. By spherical symmetry, 
the area to be minimized is 



where pm is the saddle point. In the vicinity of the horizon (y ^ 1), we can neglect the dependence 
of the functions F{y,p) and R'{y,p) on the coordinate y. The minimization of the area of the 
surface gives the equation 

= C = const. (315) 



„2 / 1 I F (dy\2 

The area of the minimal surface F/j is then given by the integral 



n2 

Area(F,,) = 2tt I dpA(p), A = , (316) 

p2 A.cv 



Using ((3TOI) we find that A{p) = + 9(2Mr+) ^ Substituting this expansion into (pT6)) we 
find that the first two terms produce divergences (when e goes to zero) which, according to our 
proposal, are to be interpreted as UV divergences of the entanglement entropy. At the black hole 
horizon, one has the relation 2Rgg\r^ — r^(i? — Raa)- Putting everything together and applying 
proposal (|304p . one finds for the divergent part 

^..= ^A^^-^/^(Ji?..4i?)lne, (317) 
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where A{J^) — Aur]^ is the horizon area. 

The logarithmic term in (|317p is is related to the logarithmic divergence (calculated holograph- 
ically in [HI]) 

in the quantum effective action of boundary CFT. This relation is a particular manifestation of the 
general formula (j280l) that relates the logarithmic term in the entropy to the conformal anomalies 
of type A and B. One notes that in the = 4 superconformal SU{N) gauge theory one has that 

It should be noted that the UV finite terms and their dependence on the size Li^v of the box 
can be computed in the limit of small Li™. This calculation is given in |201) . In particular, in any 
even dimension d one finds an universal term in the entropy that takes the form (up to numerical 
factor) S ~ r!^~^/i(d)ee(r+)i?^^, In Linv and is proportional to the value of the "obstruction tensor" 
on the black hole horizon. The direct calculation of such terms in the entanglement entropy on 
the CFT side is not yet available. 

8 Can entanglement entropy explain the Bekenstein- Hawking 
entropy of black holes? 

Entanglement entropy of a black hole is naturally proportional to the area of the black hole horizon. 
This property makes it very similar to the Bekenstein-Hawking entropy assigned to the horizon. 
This apparent similarity between the two entropies is the main motivation to raise the question, 
whether the Bekenstein-Hawking entropy is in fact entirely the entropy of entanglement. In this 
section we discuss problems which this interpretation has to face, different approaches to solve 
them and difficulties which still remain unsolved. 

8.1 Problems of interpretation of the Bekenstein-Hawking entropy as 
entanglement entropy 

Any approach that wants to treat the Bekenstein-Hawking entropy as an entanglement entropy 
has to answer the following questions: 

(i) The entanglement entropy is a UV divergent quantity, while the Bekenstein-Hawking entropy 
is a finite quantity, defined with respect to Newton's constant, measured in experiments. How can 
these two quantities be equal? 

(ii) The entanglement entropy is proportional to the number of different field species which exist 
in Nature. On the other hand, the Bekenstein-Hawking entropy does not seem to depend on any 
number of fields. This problem is known as the "species puzzle" . 

(iii) We have seen that entanglement entropy due to fields which are non-minimally coupled to 
gravity, the gauge bosons and gravitons, behave differently from the entropy due to minimally 
coupled fields. Since the gauge bosons and gravitons are fields that are clearly present in Nature 
and thus should contribute to the entropy, does this contribution spoil the possibility of interpreting 
the black hole entropy as an entanglement entropy? 

8.2 Entanglement entropy in induced gravity 

One, possibly very natural, way, originally proposed by Jacobson [142] . to attack these problems 
is to consider gravity as an induced phenomenon, in the spirit of Sakharov's ideas |189] (for a 
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review on a modern touch on these ideas see [219] ). In this approach the gravitational field is 
not fundamental but arises as a mean field approximation of the underlying quantum field theory 
of fundamental particles (the constituents). This is based on the fact, that even if there is no 
gravitational action at tree level, it will appear at one-loop. The details of this mechanism will, of 
course, depend on the concrete model. 



Model with minimally coupled fields. To start with, let us consider a simple model in which 
the constituents are minimally coupled fields: we consider Nq scalars and N1/2 Dirac fermions. 
The induced gravitational action in this model, to lowest order in curvature, is 



/ RVgd'^x, (318) 

W-KLjind J E 



where the induced Newton's constant is 

1 N 



Gind 127re2 



, = iVo + 2iVi/2 , (319) 



N is the number of field species in this model. The renormalization statement which is valid for 
the minimally coupled fields guarantees that the there is a precise balance between the induced 
Newton's constant and the entanglement entropy, so that 

= 7^^(S) = J^A{^) = Sbh , (320) 

i.e. the entanglement entropy of the constituents is precisely equal to the Bekenstein-Hawking 
entropy, expressed in terms of the induced Newton's constant pi9p . Thus, if at a fundamental 
level the constituents in Nature were only minimal fields, the Bekenstein-Hawking entropy, as this 
example shows, would be explained as the entropy of entanglement. Of course, this example ignores 
the fact that there are elementary particles, namely the gauge bosons, which are non-minimally 
coupled. 



Models with non-minimally coupled fields. In the model with minimal fields the induced 
Newton's constant (|319p is set by the UV cutoff e. If one wants to deal with the UV finite quantities 
one has to add fields which contribute negatively to Newton's constant. Excluding non-physical 
fields with wrong statistics, the only possibility is to include non-minimally coupled fields, scalars 
or vectors. Models of this type have been considered by Frolov, Fursaev and Zelnikov [M], [96] , 
[88] . One considers a multiplet of scalar fields of mass and non- minimal coupling and a set 
of massive Dirac fields with mass md- The number of fields and their parameters are fine tuned 
so that the ultra-violet divergences in the cosmological constant and in Newton's constant are 
canceled. The induced Newton's constant then 

^(l-6e.)m2lnm2 + 2^m2lnm2 j (321) 

s d ) 

is dominated by the mass of the heaviest constituents. However, as soon as we include the non- 
minimally coupled fields the precise balance between Newton's constant and the entanglement 
entropy is violated, so that the Bekenstein-Hawking entropy Sbh — ^(S)/4G'iri(ij defined with 
respect to induced Newton's constant, is no longer equal to the entanglement entropy. In the 
model considered in [94], [96], [88] (various models of a similar nature are considered in [89], [95) . 
[87] . [91], [92], [97], [101] . [107J ) the exact relation between two entropies is 

Sbh — Sent — Q , (322) 
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where the quantity Q is determined by the expectation value of the non-minimahy coupled scalar 
fields on the horizon E 

Q-2^V6 / <<t>l> . (323) 

This quantity is UV divergent. For a single field it is similar to the quantity (jl42[) . Thus the sharp 
difference between the entanglement entropy and the Bekenstein-Hawking entropy in this model 
can be summarized as follows: even though the induced Newton's constant is made UV finite, the 
entanglement entropy still (and, in fact, always) remains UV divergent. Thus, we conclude that, 
in the model of Frolov, Fursaev and Zelnikov, the entanglement entropy is clearly different from 
the Bekenstein-Hawking entropy. 



8.3 Entropy in brane-world 

An interesting example where the Bekenstein-Hawking entropy is apparently induced in the correct 
way is given in |128j . This example is closely related to the AdS/CFT correspondence discussed 
in section [71 In the Randall-Sundrum set-up [186^ one may consider the regularized boundary, 
which appeared in our discussion of section [71 as a 3-brane with Z2 symmetry in an anti-de Sitter 
space-time. In the framework of the AdS/CFT correspondence this brane has a description in 
terms of CFT coupled to gravity at a UV cutoff [125 . If the brane is placed at the distance p — 
from the Anti-dc Sitter boundary, one obtains that there is a dynamical gravity induced on the 
brane, with the induced Newton's constant 

1/Gn = 2iVV(^e') , (324) 

where TV is the number of colors in the superconformal SU{N) Yang-Mills theory. iV^ in this case 
plays the role of the number of species. We notice that, according to the AdS/CFT dictionary, the 
parameter e which is an infra-red cut-off on the anti-de Sitter side, is in fact, a UV cut-off on the 
CFT side. Consider now a black hole on the 3-brane. The Bekenstein-Hawking entropy can then 
be represented as follows 

_ A(E) _ N^AjE) 

As Hawking, Maldacena and Strominger ^128j suggested , the right hand side of p25p can be 
interpreted as the entanglement entropy of iV^ fields. This interpretation turns out to be the right 
one, if one uses the holographic entanglement entropy discussed in section [71 Indeed, taking the 
leading divergent term in (I317|) and noting that, in a Z2 brane configuration, this result should be 
multiplied by factor of 2, we get exactly the right hand side of p25p . In rcf.[128 one considers 
de Sitter spacetime (so that (I325P is the entropy of the de Sitter horizon in this case) on the 
brane since it is the simplest brane configuration one can construct in anti-de Sitter spacetime. In 
|83l this proposal was extended to the holographic entanglement entropy of black hole on 2-brane 
solution found in 84 . The two-dimensional black hole is considered in [102) . Entropy of a generic 
black hole on the 3-brane was considered in 12011. 

There are, however, certain open questions regarding this example. First of all we should 
note that the weakly coupled A/" = 4 SU{N) supermultiplet contains the Yang- Mills fields (gauge 
bosons), conformally coupled scalars and the Weyl fermions |167J . Thus it is a bit of a mystery 
how the entanglement entropy of these, mostly non-minimally coupled, fields (gauge bosons and 
scalars) has managed to become equal to the Bekenstein-Hawking entropy, when recalling the 
problems with the non- minimal coupling we have discussed in section [XTHl A part of this mystery 
is the fact that the holographic regularization (which corresponds to infra-red cut-off on the anti-de 
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Sitter side) does not have a clear analog on the boundary side. Indeed, if we take for example a 
standard heat kernel regularization, we find that the term linear in the scalar curvature R does 
not appear at all in the effective action produced by the weakly coupled J\f = 4 superconformal 
gauge multiplet. 

8.4 Gravity cut-ofF 

If we compare the two examples when the Bekenstein-Hawking entropy is correctly reproduced, 
the model of induced gravity with minimally coupled constituents and the brane world model, we 
find that the success of the two models is strongly based on a precise relation between Newton's 
constant, the number of species and the UV cut-ofF. This relation can be reformulated in terms of 
the Planck mass MpL (1/Gjv ^ ^J^pl) ^^"^ energy cutoff A '-^ 1/e, 

A = Mpl/Vn, (326) 

where the precise numerical pre-factor depends on the concrete model. It is amazing to note that 
exactly this relation (|326|) was proposed by Dvali [77], [79], [78] to hold in general in a theory of 
Quantum Gravity coupled to a large number of matter species. The arguments which were used 
to get this relation, although they include some thought experiments with black holes, are, in fact, 
unrelated to (and are thus independent of) the entropy. However, it is clear that relation p26|) . 
provided it is correct, helps to reproduce precisely the Bekenstein-Hawking entropy as entropy of 
entanglement and automatically solves the species puzzle [80] , 

8.5 Kaluza-Klein example 

One example when relation (I326P holds is the Kaluza-Klein model. In this model one starts with 
(4 -t- n)-dimensional theory of gravity which then is compactified so that one has n compact direc- 
tions, forming, for example, an n-tours with characteristic size R, and 4 non-compact directions 
which form our 4-dimensional geometry. The higher-dimensional Planck scale A is considered fun- 
damental in this model and plays the role of the UV cutoff while the 4-dimensional Planck scale 
MpL (or 4-dimensional Newton's constant) is derived, 

Af|,i = A2(i?A)" . (327) 

Suppose that in higher dimensions there is only one particle - the massless graviton. From the 
four-dimensional point of view one has, additionally to a single massless graviton, a theory of the 
tower of spin-2 massive Kaluza-Klein (KK) particles. Truncating the tower at the cut-off A one 
finds that N = (_RA)" is precisely the number of these Kaluza-Klein species. Thus, as was noted 
in [77], [79], the relation (|327p is a particular example of relation (|326p in which N should be 
understood as the number of the KK species. 

In the Kaluza-Klein example the entanglement entropy is equal to the Bekenstein-Hawking 
entropy as demonstrate in Dvali and Solodukhin |80j . Consider now a large black hole with horizon 
size rg ^ R. This large black hole fills all compact directions so that, from the higher-dimensional 
point of view, the black hole horizon is a product of a 2-sphere of radius Vg and an n-dimensional 
torus of size R. The Bekenstein-Hawking entropy in the (4 -I- n)-dimensional theory is 

S-jj^+'O ^ 4^A"+2r2i?" , (328) 

where AirrgK" is the area of (4 -t- n)-dimensional horizon. From the 4-dimensional point of view 
this horizon is a 2-sphere of radius rg and the Bekenstein-Hawking entropy in the 4-dimensional 
theory is 

S'j^l = 4nMj,^rl . (329) 
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We remark observe that these two entropies are equal so that the two pictures, the higher dimen- 
sional and 4-dimensional one, are consistent. Let us discuss now the entanglement entropy. In the 
(4 + n)-dimensional theory there is only one field, the massless graviton. Its entropy is 

^ 4^r2^«A"+2 , (330) 

where the cut-off is the higher dimensional Planck scale A. On the other hand, in the 4-dimensional 
theory one computes the entanglement entropy of N KK fields 

Si^j, = NiinrlA') = AnrlA' (i?A) . (331) 

These two entropies are equal to each other so that the two ways to compute the entangle- 
ment entropy agree. Moreover the entanglement entropy (|33ip and (|330|) exactly reproduces the 
Bekenstein-Hawking entropy p28|) and p29|) . 

Discussing this result we should, however, note that the massless and massive gravitons are non- 
minimally coupled particles. It remains to be understood how the problem of the non-minimal 
coupling is overcome in this example. 

9 Other directions of research 

In this section we briefiy mention some other interesting directions of on-going research. 
9.1 Entanglement entropy in string theory 

It is generally believed that the problem of entanglement entropy of black hole can and should 
be resolved in string theory. This was originally suggested by Susskind and Uglum [212) . Indeed, 
assuming that string theory is UV finite, the corresponding entropy calculation should result in 
a finite quantity. More specifically, the effective action of a closed string can be decomposed 
in powers of string coupling g as 5^'""^', where n is the genus of the string world sheet. The 
string configurations with spherical topology, n = 0, give a contribution. In a low energy 
approximation this is exactly the contribution to Newton's constant G ^ . Thus one may expect 
that taking into account just n = closed string configurations, one will correctly reproduce both 
the entanglement entropy and Newton's constant. In the Euclidean formulation the prescription 
of [212] is to look at the zero genus string world sheet which intersects the Killing horizon. In 
the Lorentzian picture this corresponds to an open string with both ends attached to the horizon. 
The higher genus configurations should give some corrections to the n = result. This is a very 
attractive idea. However a very little progress has been made in the literature to actually calculate 
the entanglement entropy directly in string theory. The reason is of course the technical complexity 
of the problem. Some support to the idea of Sussking and Uglum was found in work of Kabat, 
Shenker and Strassler |151j . where the entropy in a two-dimensional 0{N) invariant cr-model was 
considered. In particular, it was found that the state counting of the entropy in the UV regime 
may be lost if considered in the low energy (IR) regime. This type of behavior models the situation 
with the classical Bekenstein-Hawking entropy. Presumably this analysis could be generalized to 
the string theory u- model considered either in optical target metric [9], [10] or in the Euclidean 
metric with a conical singularity at horizon (as suggested in [33]). Possibly in the latter case the 
results obtained for strings on orbifolds ^ can be useful (see [50], [IH], [SJ for earlier works in 
this direction). 

Another promising approach to attack the problem is to use some indirect methods based 
on dualities. For example, the AdS/CFT correspondence has been used in p5] to relate the 
entanglement entropy of a string propagating on gravitational AdS background with a Killing 
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horizon to the thermal entropy of field theory defined on a boundary of AdS and then, eventually, 
the thermal entropy to the Bekenstein-Hawking entropy of the horizon. 

An interesting approach to the entanglement entropy of extremal black holes via AdS2/CFTi 
duality is considered in [4], [191) . where, in particular, one can identify the entanglement entropy 
and the microcanonical statistical entropy. This approach is based on the earlier work of Maldacena 
|168] in which the Hartle-Hawking state is identified with an entangled state of two copies of CFT, 
defined on two boundaries of the maximally extended BTZ spacetime. In the accurately taken 
zero temperature limit the reduced density matrix, obtained by tracing over the states of one copy 
of CFT, of the extremal black hole is shown to take the form 

d{N) 

^=d(]^El^><^l' (332) 

which describes the maximally entangled state in the two copies of the CFTi living on the two 
boundaries of global AdS2 . d{N) is the dimension of the Hilbert space of CFTi , it can be expressed 
in terms of the charges carried by the black hole. The corresponding entanglement entropy S = 
— Trplnp = In (i(iV) then is precisely equal to the microcanonical entropy in the familiar counting 
of BPS states and thus is equal to the black hole entropy [32] . 



9.2 Entanglement entropy in loop quantum gravity 

Another approach to Quantum Gravity, sometimes considered as competing with string theory, is 
Loop Quantum Gravity. In this theory one considers polymeric excitations of the gravitational 
field represented by the states of spin networks. A spin network is a graph, a network of points 
with links representing the relation between points. Each link is labeled by a half-integer j (the 
label stands for SU{2) representations). To points, or vertices, of a spin network are attached a 
SU{2) so-called intertwiner, a SU{2) invariant tensor between the representations attached to all 
the edges linked to the considered vertex. A simpler and more familiar object in particle physics is 
the Wilson loop. A surface S is represented by vertices (punctures) which divide the spin network 
on two parts. By tracing over states of just one part of the network one obtains a density matrix. 
The entanglement entropy then reduces to a sum over intersections of the spin network with the 
surface E [56], [165] . [66] . 

p 

5(S]) = ^ln(2jp + l), (333) 

where P is number of punctures representing E. This quantity should be compared to eigenvalues 
of the operator of area, 

A(E) = 8^G7^V^jp(jp + l). (334) 
p=i 

Both quantities scale as P for large P, which indicates that the area law is correctly reproduced. 
The exact relation between the two quantities and the classical entropy Sbh — ^(E)/G is, however, 
not obvious due to ambiguities present in the formalism. The Immirzi parameter 7 is one of them. 
The question, whether the Bekenstein-Hawking entropy is correctly reproduced in this approach, 
is eventually related to the continuum limit of the theory |143j . As discussed by Jacobson 143], 
the answering this question may require certain renormalization of Newton's constant as well as 
area renormalization. Indeed, quantity p34p represents a microscopic area which may be related 
to the macroscopic quantity in a non-trivial way. These issues remain open. 
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9.3 Entropy in non-commutative theories and in models with minimal 
length 

One might have hoped that the UV divergence of the entanglement entropy could be cured in a 
natural way were the structure of spacetime modified on some fundamental level. For example, if 
spacetime becomes non-commutative at short distances. This idea was tested in the case of simple 
fuzzy spaces in [68] . [67]. Although the area law has been verified, the entanglement entropy 
appears to be sensitive to the size of the ignored region, a phenomenon which may be understood 
as a UV-IR mixing typical for the non-commutative models. 

A holographic calculation of the entanglement entropy in non-commutative Yang-Mills theory 
was considered in [M], [13]. This calculation for a strip of width / shows that for for large values 
of I ^ Ic compared to some characteristic length Ic ~ flA^/^/e, where is the parameter of non- 
commutativity and A = gyM-^ the 't Hooft coupling, then the short-distance contribution to 
the entanglement entropy shows an area law of the form 

S-N^ff^, N.ff=NH^-^), (335) 

while for smaller values I ~ Ic the entropy occurs to be proportional to the volume. As seen 
from p35p the non-commutativity does not improve the UV behavior of the entropy but leads to 
the renormalization of the effective number of degrees of freedom that may be interpreted as a 
manifestation of non-locality of the model. 

The other related idea is to consider models in which the Heisenberg uncertainty relation is 
modified as AxAp > |(1 -I- X'^{Ap)^), which shows that there exists a minimal length Aa; > hX 
(for a review on the models of this type see [113] ). In a brick wall calculation the presence of this 
minimal length will regularize the entropy as discussed in [27], [224] . |209j . |155) . [156] . 

9.4 Transplanckian physics and entanglement entropy 

One way to check whether the entanglement entropy is sensitive to the way the conventional theory 
is completed in the UV regime is to study the possible modifications of the standard Lorentz 
invariant dispersion relation = at large values of momentum k (or at short distances). A 
typical modification is to break the Lorentz invariance as follows oj^ = k-^ + /(k^). This issue was 
studied in [147J and f40] in the context of the brick wall model. The conclusions made in these 
papers are however opposite. According to ref.^l47j the entropy is still UV divergent although 
the degree of divergence is modified in way which depends on the form of function /(k^). On the 
other hand, the paper [ID] claims that the entropy can be made completely UV finite. In a similar 
claim ref. [185J suggests that the short-distance finiteness of the 2-point correlation function should 
imply the UV finiteness of the entanglement entropy. The entanglement entropy in a wide class 
of theories characterized by modified (Lorentz invariant or not) field operators (so that the UV 
behavior of the modified propagator is improved compared to the standard one) was calculated 
in |183| . The conclusion reached in [183] (see also discussion in sections 12.121 and 13.131 of this 
review) agrees with that of ref. |147) : no matter how well is the UV behavior of the propagator, 
the entanglement entropy remains UV divergent. That the short-distance regularity of correlation 
functions does not necessarily imply that the entanglement entropy is UV finite was pointed out 
in [182] . 

9.5 Entropy of more general states 

The quantum pure state which is the starting point in the entanglement entropy calculation should 
not necessarily be a vacuum state. According to refs.[SS], [SI], [S5], [S3], if one starts with a mixed 
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state of the vacuum and an excited state, the entanglement entropy receives some power law 
corrections, 

5^4(l + c(4)'^)' (336) 
whith /? always less than unity, and the power law correction is due to the excited state. 

9.6 Non-unitary time evolution 

An interesting issue discussed in the literature is the time evolution of entanglement entropy. It 
was suggested in [2^ that the eigenvalues of a reduced density matrix depend on time t. This is 
not possible if the time evolution of the density matrix is described by an unitary operator. Thus 
the time evolution should be nonunitary. In particular, the entanglement entropy should depend 
on time t. Similar conclusions have been made in [31] and |158) . jl59) . [160J . where, in particular, 
it was shown that the entanglement entropy is an increasing function of time. These observations 
may have interesting applications for black holes. As was proposed by Hawking [130] the evolution 
in time of a black hole should be nonunitary, so that a pure initial state may evolve into a mixed 
state. From the entanglement point of view, this behavior appears to be not in contradiction with 
the principles of quantum mechanics, rather it is a simple consequence of the entangled nature of 
the system. The irreversible loss of information due to entanglement is also seen from the evolution 
of the entropy under RG flow ^U^, \TM\ . 

10 Concluding remarks 

Since the inspiring work of Srednicki in 1993 we have come a long way in understanding the 
entanglement entropy. Of course, the main motivation for this research always was and still is the 
attempt to consider the Bekenstein-Hawking entropy as entropy of entanglement. If successful this 
approach would give a universal explication for the entropy of black holes, valid for black holes of 
any size and mass and carrying any charges. Unfortunately, we are not yet at that point. Many 
important unresolved problems remain. However, after 17 years of continuous progress it would 
not be too surprising if we were actually not that far from the final answer. Speaking about the 
future developments, I think that the further progress could be made in advancing in two main 
directions: understanding the entanglement entropy directly in string theory and resolving the 
puzzle of the non-minimal coupling. Hopefully, at some not so distant point in time the efforts 
made in these directions will provide us with the still missing elements in the picture. 
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